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ABSTRACT 


Optimum  reduction  techniques  have  been  determined  and  a  computer 
;  program  written  to  measure  density,  winds,  temperature  and  pressure  from 

the  passive  ROBIN  Sphere  in  the  30-100  km  region  of  atmosphere.  The 
program  extends  the  altitude  range  of  the  1965  ROBIN  Program  from  70  to 
'  100  km.  Output  from  the  computer  program  also  includes  the  estimated 

noise  errors  in  density,  winds,  temperature,  and  pressure,  the  frequency 
response  of  density  and  winds  as  a  function  of  altitude,  and  Cal -comp  plots 
of  the  meteorological  parameters. 

The  smoothing  techniques  used  in  the  program  were  determined  so 
as  to  minimize  the  sum  of  the  random  and  bias  errors  in  density  and  winds. 
A  separate  smoothing  function  was  determined  for  density  and  for  wind 
calculations.  The  density  smoothing  function  expands  below  60  km  in  order 
to  maintain  a  constant  2%  noise  error  in  density. 

All  possible  sources  of  error  in  the  meteorological  parameters  are 
discussed.  For  the  first  time  error  equations  have  been  derived  for 
temperature  and  pressure  errors.  Error  estimates  are  made  for  random 
errors,  biar  errors,  and  errors  resulting  from  the  initial  estimate  of 
temperature. 
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1. 


INTRODUCTION 


1.1  Areas  ROBIN  System 

The  ROBIN  system  consists  of  a  ROBIN  sphere,  an  Areas  rocket 
motor  and  an  AN/FPS-I6  tracking  radar.  The  ROBIN  sphere  is  made  of 
one  half  mil  Mylar  inflatable  to  a  diameter  of  one  meter  containing  an 
internally-supported  corner  reflector.  Packaged  in  a  collapsed  condition 
within  the  nosecone  of  a  meteorological  rocket,  it  is  ejected  at  the  apogee  of 
the  rocket  and  inflated  by  vaporization  of  isopentane,  to  a  super  pressure 
of  approximately  ten  millibars.  Thus  inflated,  the  ROBIN  sphere  is  tracked 
from  apogee  to  approximately  30  kilometers  by  an  AN/FPS-16  high  precision 
tracking  radar.  The  Areas  rocket  motor  is  a  4.  5  inch  diameter  solid  pro¬ 
pellant,  end  burning  rocket  capable  of  carrying  the  sphere  payload  to  an 
altitude  of  75  kilometers.  The  FPS-16  tracking  radar  generates  spherical 
space -time  coordinates  at  digitized  increments  of  1/10  of  a  second.  From 
the  space-time  coordinates,  the  meteorological  parameters  of  density,  wind, 
temperature  and  pressure  are  deduced.  A  complete  discussion  of  the  Areas 
ROBIN  system,  its  advantages  and  shortcomings  is  contained  in  the  report 
by  Engier,  "Development  of  Methods  to  Determine  Winds,  Density,  Pressure 
and  Temperature  from  the  ROBIN  Falling  Balloon,  1965".  (Ref.  1). 

1.  2  Viper  Dart  ROBIN  System 

In  recent  years  the  Air  Force  has  extended  the  concept  of  the  falli  ig 
sphere  experiment  to  measurements  at  altitudes  from  70  to  100  km.  To 
implement  this  concept,  the  capability  of  the  rocket  vehicle  had  to  be  increased 
to  permit  sphere  apogees  in  excess  of  120  km.  Boosted-Arcas  vehicles  have  been 
utilized  but,  because  of  their  high  cost,  have  generally  been  superceded  by  a  lower 
cost  vehicle,  the  Viper  Dart,  developed  especially  for  the  small  light-weight 
ROBIN  payload.  The  Viper  Dart  rocket  (Ref.  II)  was  developed  successfully  and 
has  been  employed  to  carry  the  sphere  to  its  required  apogee.  Although  the 
March  1965  ROBIN  Computer  Program  gave  highly  accurate  measurements  below 
70  km,  it  was  not  satisfactory  for  measurement  above  70  km.  (see  Ref.  2).  Hence 
a  new  high  altitude  ROBIN  program  was  required  for  measurements  extending  to 
100  km.  The  problems  encountered  with  the  1965  program  which  had  to  be  dealt 
with  in  designing  a  high  altitude  program  were: 

i)  The  smoothing  used  in  the  March  1965  ROBIN  program  introduced 
a  significant  bias  in  density  above  70  km. 

ii)  Above  80  km,  wind  accuracy  decreased  rapidly  with  increasing 
altitude. 

iii)  An  examination  of  the  adequacy  of  the  Lambda  (X )  check  below  40  km 
was  in  order.  A  high  altitude  collapse  check  also  had  to  be  deter¬ 
mined. 
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It  is  the  purpose  of  this  report  to  discuss  the  new  high  altitude  data 
reduction  program  which  is  titled,  "The  May-1970  AFCRL  ROBIN  Program", 
explain  the  rationale  and  methodology  used  to  design  the  program,  and  to  discuss 
the  errors  in  the  winds  and  thermodynamic  data  that  results  with  the  use  of  this 
program, 

l.  2. 1  Program  Specifications 

The  preliminary  specifications  for  the  May-1970  AFCRL  ROBIN  Program 
consisted  of  the  following  items:  a)  the  program  should  be  optimum  for  measuring 
density  and  wind  in  the  70  to  100  kilometer  region  of  the  atmosphere,  assuming  a 
sphere  apogee  of  125  kilometers;  b)  the  program  should  also  give  accurate  and 
reliable  density  measurements  from  30  to  70  kilometers;  c)  even  though  the  data 
reduction  technique  need  not  be  optimum  for  sphere  apogees  other  than  125  kilo¬ 
meters,  other  sphere  apogees  between  75  and  140  kilometers  should  not  result 
in  a  serious  degradation  of  the  accuracy  of  the  meteorological  parameters; 
d)  temperature  and  pressure  accuracies  should  be  commensurate  with  density 
accuracy;  e)  the  program  should  accurately  determine  the  altitude  of  sphere 
collapse. 

1.  2.  2  Discussion  of  Coordinate  System 

The  coordinate  system  chosen  is  a  left-handed  system  in  which  the 
two  horizontal  directions  are  X  (down-range)  and  Y  (cross-range).  Looking 
down -range,  the  positive  cross -range  direction  is  to  the  right.  The  positive 
clockwise  angle  from  geographical  north  to  the  positive  X-axis  of  the  range 
coordinate  system  is  called  Z  (see  Figure  1).  The  azimuth  angle  (A)  is  the 
angle  clockwise  from  the  positive  X-axis.  The  vertical  direction  is  Z,  and 
its  negative  direction  is  toward  the  center  of  the  earth  at  the  origin  of  the 
coordinate  system.  Hence,  the  horizontal  X-Y  plane  (assuming  a  spherical 
earth)  is  tangent  to  the  earth  at  the  launch  site.  The  altitude  (h)  is  defined 
as  the  perpendicular  distance  from  the  sphere  to  the  earth. 


Figure  1:  Three-dimensional  Coordinate  System  used  in  Defining  Equations 
oi  Motion. 
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In  contrast  to  the  coordinate  system  defined  above,  the  meteorologist 
desires  wind  velocities  with  respect  to  the  coordinate  system  with  "h"  defined 
as  the  vertical  axis.  To  determine  the  practical  effect  of  the  differences 
between  the  two  coordinate  systems,  winds  were  computed  from  several 
flights  using  each  coordinate  system.  The  largest  difference  in  component 
wind  velocity  between  the  two  systems  was  0.  2  m/sec.  Winds  are  initially 
computed  in  the  original  coordinate  system  (with  Z  as  vertical  axis)  and  later 
transformed  to  the  coordinate  system  with  Mh"  defined  as  the  vertical  axis. 

1.2.  3  Density  and  Wind  Measurement 

To  obtain  density,  the  drag  force  that  the  atmosphere  exerts  upon 
the  sphere  must  be  measured.  In  the  70  to  100  kilometer  region,  the  vertical 
velocities  and  accelerations  are  much  larger  than  the  horizontal  velocities 
and  accelerations.  For  this  reason  the  drag  acceleration  is  primarily  in  the 
vertical  direction.  Accurate  density  calculations  are  thus  largely  a  result  of 
the  accuracy  to  which  vertical  velocities  and  accelerations  can  be  measured. 

Horizontal  winds  influence  the  sphere's  trajectory  by  inducing 
horizontal  excursions  in  its  path  in  three  dimensional  space.  These  horizontal 
excursions  are  used  to  reconstruct  the  wind  profile.  At  altitudes  above  85  km, 
the  vertical  acceleration  also  plays  an  important  role  in  reconstructing  the 
horizontal  wind  profile.  Thus,  for  measuring  wind,  the  horizontal  velocity 
and  acceleration  components,  as  well  as  the  vertical  acceleration,  must  be 
determined  accurately. 

Since  density  is  derived  primarily  upon  vertical  measurements  and 
winds  depend  on  vertical  as  well  as  horizontal  measurements,  two  separate 
smoothing  functions  have  been  introduced  into  the  program.  The  first 
smoothing  function  is  defined  so  as  to  optimize  density  measurements.  This 
smoothing  is  applied  only  to  vertical  position  coordinates  to  obtain  vertical 
velocities  and  accelerations.  A  second  smoothing  function  is  defined  so  as  to 
optimize  wind  measurements.  This  smoothing  is  applied  to  horizontal 
position  coordinates  to  obtain  horizontal  velocities  and  accelerations  as  well 
as  reapplying  it  to  vertical  position  to  obtain  a  second  estimate  of  vertical 
velocity  and  acceleration  to  be  used  in  the  wind  equation.  The  determination 
of  these  two  smoothing  functions  is  the  critical  problem  in  devising  a 
reduction  program.  The  rationale  and  methodology  used  to  define  these 
smoothing  functions  is  set  forth  in  the  following  sections. 

2.  DENSITY 

Density,  p,  is  computed  by 


m(g  -  Z  -  C  )  (1) 

z  z 

1/2  C  A  V(Z  -W  )  +  VLg  (See  Ref.  1) 

D  z  b  z 
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where 


Mass  of  sphere 

Gravitational  acceleration  in  vertical  direction 

Vertical  acceleration  of  sphere 
Coriolis  acceleration  in  vertical  direction 

Drag  coefficient 

Cross  seciion  area  of  sphere 
Velocity  o;!  the  sphere  relative  to  the  air 
Vertical  velocity  of  sphere 
Vertical  wind  velocity 

Volume  of  sphere 


The  computed  density  error  is  a  result  of  the  errors  present  in  the 
parameters  c*n  ihe  right  side  of  Equation  1.  The  error  present  in  many  of  the 
variables  of  Equation  1  makes  a  negligible  contribution  to  the  error  in  density. 


2.  1  Variables  Not  Producing  Significant  Density  Errors 


The  vertical  component  of  the  acceleration  due  to  gravity  (g  )  is 
computed  by  the  equation 


gs  (Z+r) 
(l+j;)2  (h+r) 


(2) 


« 


► 


* 


where 

g  =  Sea  level  gravitational  constant 

r  =  Radius  of  the  earth 

h  =  Altitude  of  sphere 

Z  =  Vertical  position  of  sphere 

The  accuracy  of  Equation  2  is  sensitive  only  to  the  accuracy  of  the 
radius  of  the  earth,  or  more  precisely,  the  determination  of  the  center  of 
mass  of  earth.  Even  if  the  radius  to  the  center  of  mass  of  the  earth  were 
in  error  by  10  kilometers,  the  resultant  error  in  density  would  be  less  than 
1/2%  at  an  altitude  of  100  kilometers  and  much  less  at  lower  altitudes. 

The  mass  of  the  sphere  (m)  is  determined  for  each  sphere  individually  * 

by  the  manufacturer.  The  accuracy  of  the  mass  thus  obtained  should  be  to  at 
least  1/2  gram.  For  a  nominal  spnere  mass  of  120  grams  this  would  result 
in  less  than  1/2%  density  error. 
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Manufacturer's  specifications  for  the  ROBIN  sphere  demand  that  it  should 
deviate  less  than  1%  from  a  perfect  sphere.  A  1%  error  in  diameter  produces  a  2% 
error  in  A  and  thus  a  2%  bias  in  density.  An  error  in  either  diameter  or  mass  would 
not  effect  the  shape  of  the  density  profile  but  only  translate  it  in  some  direction.  This 
results  because  the  error  in  diameter  or  mass  is  constant  for  the  entire  flight.  Thus, 
to  a  first  order  approximation,  a  constant  bias  results.  However,  if  mass  should  vary 
due  to  a  gas  leak,  the  bias  error  would  not  be  constant. 

From  30  to  100  kilometers  the  bouyancy  (V^g  )  and  coriolis  (C  )  terms 
make  only  minor  contributions  to  the  determination  of  density.  Hence,  any  error 
present  in  these  terms  gives  negligible  contribution  to  density  error. 

The  remaining  variables  which  could  make  a  significant  contribution  to 

density  error  are  C  .  4,  2,  W  ,  and  V. 

D  z 


2.  2  Density  Error  Equation 


For  the  purpose  of  deriving  an  error  equation  for  density,  the  density 
equation  can  be  simplified  to  Equation  3 


m{g  -Z) 
z 

1/2  C  A  V  (Z-W  ) 

D  Z 


(3) 


where  both  the  bouyancy  and  coriolis  forces  have  been  neglected. 


Considering  the  error  in  density  to  be  a  function  of  the  errors  in 
Cjy  Wz,  Z,  Z,  and  V.  The  error  increment  in  density  due  to  errors  in  the  five 

parameters  is  given  by  the  differential  approximation  formula  as 


Ap  = 


dp 


bC 


AC 


D 


bP. 
D  '  dV 


AV  +  -S- 
d  Z 


AZ  + 


dp 

dlT 


AZ  + 


bP 

'dW 


AW 


(4) 


Each  of  these  six  terms  can  be  analyzed  to  determine  the  magnitude  of  its 

contribution  to  the  error  in  density.  Since  Cj-j,  Wz  and  errors  in  V  (primarily  a 

result  of  horizontal  wind  errors)  are  essentially  independent,  they  can  be  analyzed 

•  •  • 

individually.  Errors  in  Z  and  Z  are  not  independent  due  to  the  smoothing  process. 
Thus,  the  terms  containing  2  and  Z  must  be  considered  together.  The  following 
sections  describe  the  error  terms  and  provide  means  of  estimating  their  effect  on 
density  accuracy. 


2.  2. 1  Drag  Coefficient  Error  Term 


The  percent  error  in  density  in  terms  of  the  percent  error  in  C  is 

given  by  the  error  variance  formula  tr  ^  ^  ^ 

P 


C 

^  —  y  (See  Appendix  A) 
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а)  Error  due  to  Drag  Table 

Since  the  accuracy  of  the  drag  coefficients  cannot  be  determined 
precisely,  it  is  impossible  to  give  an  exact  error  variance  value  for 
the  resulting  per  cent  error  in  density.  The  accuracy  of  any 
individual  drag  table  in  the  literature  today  is  suspect.  There  are 
several  experimenters,  Aroesty  {Ref.  3),  Ashkanei  (Ref.  4), 

Heinrich  et.  al.  (Ref.  5  and  6),  and  Goin  and  Lawrence  (AEDC) 

(Ref.  7)  who  have  measured  drag  coefficients  for  a  passive  sphere. 

The  inconsistencies  of  results  obtained  by  these  experimenters  have 
been  pointed  out  by  Luers  (Ref.  8)  and  others.  Until  further 
concentrated  effort  is  supported  by  interested  parties,  the  accuracy 
cf  any  drag  coefficient  presently  in  use  is  not  known,  and  probably 
varies  for  different  sections  of  the  tables.  As  a  result,  the  density 
error  due  to  drag  coefficient  error  cannot  be  determined  directly 
and  will  probably  vary  over  different  segments  of  the  flight. 

The  drag  table  to  be  used  in  the  High  Altitude  ROBIN  Program 
appears  as  figure  2  and  is  basically  the  work  of  Heinrich  (Ref.  5  and 

б)  in  the  supersonic  section,  and  of  Arnold  Engineering  Develop¬ 
ment  Center  (Ref.  7)  in  the  subsonic  region.  The  low  Reynolds 
number  (<50)  high  Mach  number  drag  coefficients  were  taken  from 
Potter  and  Miller  (Ref.  9).  The  impressive  aspect  of  this  drag  table 
is  the  similar  shapes  of  the  C^  curves  given  as  a  function  of  Mach 
and  Reynolds  number  even  though  the  drag  table  was  the  result  of  two 
independent  researchers  using  two  different  techniques  for  calculating 
drag.  (Heinrich  used  a  stationary  sphere  in  a  wind  tunnel;  Lawrence 
and  Goin  used  a  ballistic  range).  However,  even  though  this  drag 
table  shows  smooth,  consistent  drag  curves,  it  is  impossible  to  quote 
specific  accuracies  of  the  drag  table  because  of  the  interpolated  section 
of  the  table.  The  stated  accuracies  by  the  experimenters  are  as 
follows:  Lawrence  and  Goin  subsonic,  approximately  2%;  Heinrich 
supersonic;  maximum  possible  error  ranges  from  ±2.  37%  to  ±27.  9% 
depending  on  Reynolds  number.  However,  actual  errors  are  usually 
not  the  maximum  possible.  By  interpreting  Heinrich's  maximum 
possible  error  as  a  3cr  error  and  by  extrapolating  Heinrich's  and  AEDC 
error  values  in  the  interpolated  section  of  the  drag  table  a  crude 
approximation  to  Cq  error  as  a  function  of  sphere  altitude  was 
derived  for  a  sphere  with  apogee  of  125  km.  The  results  are  as 
presented  in  Table  1. 

As  of  this  writing  a  very  comprehensive  sphere  drag  determination 
program  is  being  undertaken  at  AEDC  for  AFCRL  which  will  cover 
broad  areas  of  interest  for  both  the  ROBIN  and  accelerometer  equipped 
spheres.  When  the  tests  are  completed,  it  is  expected  that  further 
analysis  will  indicate  the  need  to  modify  the  values  of  Cjj  used  herein. 


6 


_J - -J - *- 

05  10  15 


DRAG  TABLE 

SUBSONIC  -  AEDC,  LAWRENCE  and  GOIN 
SUPER  SONIC  -  HEINRICH 
LOW  REYNOLDS  -  AEDC,  POTTER  and  MILLER 


2.0  25  3.0 

MACH  NUMBER 


_j - 1 - L~ 

35  40  45 


Figure  2:  Drag  Curves  used  in  High  Altitude  ROBIN  Program. 
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Figure  2:(Continued) 
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If  so  a  new  table  of  C  will  be  incorporated  into  the  High  Altitude 
Program.  At  that  time  the  error  in  density  due  to  the  error  in 
will  be  re-evaluated. 

table  I 

Approximation  to  Percent  Error  in  Cp  -vs  -  Attitude 
for  a  Sphere  with  Apogee  12S  Km 


Sphere  Altitude  "“Km  %  Error  in  Cn 


100 

95 

90 


10 

8 

6 


85 


3 


80 


1 


75  2 

70  6 

65  4 

60  3 


55  and  below 


2 


b)  Error  due  to  Initial  Estimate  of  Temperature 


In  addition  to  the  error  in  density  due  to  the  inaccuracies  of  the 
drag  table,  a  density  error  also  results  from  the  error  in  the  initial 
estimate  of  temperature. 

At  the  first  point  of  density  computation  an  initial  estimate  of 
temperature  (or  pressure)  is  required  (see  Ref.  1).  The  error  in  the 
initial  temperature  produces  errors  in  succeeding  temperatures  that 
decrease  in  magnitude.  These  temperature  errors  produce  errors 
in  Mach  number  and  Reynold’s  number  since  both  are  functions  of 
temperature.  Consequently,  an  error  in  and  density  results,  xhe 
magnitude  of  the  resultant  density  error  depends  upon  the  magnitude  of 
.ie  initial  error  in  temperature  and  the  shape  of  the  drag  curves.  The 
shape  of  the  drag  curves  is  important  since  this  determines  the  change 
in  resulting  from  a  given  error  in  Mach  and  Reynold's  numbers. 


The  density  error  profile  due  to  a  10%  error  in  the  initial 
estimate  of  temperature  is  presented  in  Table  2.  Table  2  was 
derived  by  first  determining  the  trajectory  of  the  ROBIN  Sphere 
assuming  it  fell  in  the  1962  Standard  Atmosphere  (i.  e. ,  Standard 
Atmosphere  density  (p,,^,,)  and  temperature  (T,,^,,)  from  a  125  km 
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apogee.  Using  this  trajectory  as  if  it  represented  a  radar  track, 
temperature,  density  and  pressure  were  computed  using  the  1965 
ROBIN  Program.  A  ±10%  error  in  the  initial  estimate  of  temperature 
was  introducrtd  by  taking  Tq  =  1.  IT»«£2"  *nd  Tq  =  .  9Tn^2"  a*  *n***** 
point  of  computation.  As  seen  from  Table  2,  the  resultant  error  in 
density  becomes  less  than  1%  after  approximately  4  km  of  flight. 


Table  2: 

Percent  Errc*  in  Density 

Resulting  from  s  10%  Error  in  the  Initial  Temperature 
ALTITUDE  -  KM 


100 

95 

90 

85 

<85 

L<loq 

4.0 

.8 

.3 

.  1 

Less  than  0. 1% 

2.  2.  2  Relative  Velocity  (V)  Error  Terms 


The  percent  error  variance  of  density  as  a  function  of  the  error  in 

.2  x  /v  wy) 2  +  (Z  -  Wz)Z  1  1/2 


relative  velocity  y=  [  (X  -  Wxy  +  (Y 


is 


( 


( 


X-W 

x 


■) 


(5) 

(See  Appendix  A) 


e  • 

The  variables  X,  Y,  W  ,W  and  V  are  basically  dependent  upon 

a)  the  sphere  apogee,  b)  the  horfzon^al  velocity  imparted  to  the  sphere  by  the 

rocket,  and  c)  the  wind  field  experienced  by  the  sphere.  Since  these  variables 

change  from  flight  to  flight.  Equation  5  was  evaluated  using  the  values  of  these 

variables  from  a  representative  flight  with  sphere  apogee  125  km.  <r  and  <r 

w  w 

were  each  taken  to  be  10  m/sec  above  70  krn  and  less  at  lower  altitudes.  Late/ 
in  this  report  it  will  be  shown  that  the  mean  wind’  profile  can  be  measured  to  an 
accuracy  in  excess  of  that  used  in  deriving  Table  3.*  Table  3  shows  the  percent 
error  in  p  resulting  from  horizontal  wind  errors.  The  table  indicates  that  the 


4 

* 


♦If  a  high  frequencyfwavelength  <5  km)  large  amplitude  (in  excess  of  10  m/seef 
wind  structure  exists  above  70  km,  then  the  smoothing  filter  used  for  determining 
winds  will  not  see  these  waves  and  actual  wind  errors  may  exceed  10  m/sec 
for  these  high  frequency  waves.  However,  such  large  wind  errors  will  not  result 
in  the  density  errors  that  would  be  predicted  by  Equation  5,  since  the  density 
filter  will  compensate  by  not  allowing  the  resulting  ficticious  density  waves 
derived  from  the  wind  error  to  pass. 
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maximum  error  in  density  is  only  1. 1%  and  occurs  at  70km.  Actually,  the 
1. 1%  estimate  at  70  km  is  undoubtedly  too  large  since  the  horisontal  wind  errors 
at  70  km  are  generally  much  less  than  10  m/sec.  Hence  it  can  safely  be 
stated  that  for  the  entire  altitude  range  from  30  -  100  km,  the  error  in  density 
resulting  from  the  inaccuracies  of  the  horizontal  wind  measurements  is  less 
than  1%. 

Equation  5  also  shows  the  need  for  computing  horizontal  winds  at  all 
altitudes  at  which  density  is  computed  (i.  e.  Density  calculations  should 
commence  at  the  serae  altitude  as  horizontal  winds  or  higher).  For  if  hori¬ 
zontal  winds  are  not  computed  and  taken  to  be  zero,  a  100  m/sec  horizontal 
wind  could  induce  a  10%  error  in  densitv. 


TABLE  3 

Density  Error  Resulting  from  Horizontal  Wind  Error*  of  10,  5, end  2  m/*e c 


Altitude 


from  VIPER  DART  13 


%  Density 
Wind  Error  I  Error 


(op  If  )(100) 


Z.  Z.  3  Vertical  Velocity  and  Acceleration  Error  Terms 


Two  types  of  errors  are  present  in  velocity  and  acceleration 
measurements  when  a  smoothing  process  is  employed  to  obtain  these 
measurements:  noise  error  and  bias  error.  The  noise  error  is  the  error 
resulting  from  the  noise  present  in  the  radar  coordinates.  A  bias  error 
results  when  the  smoothing  function  does  not  adequately  represent  the  true 
path  of  the  sphere.  The  error  in  density  due  to  the  noise  and  bias  errors  in 
velocity  and  acceleration  is  given  by  the  terms 


( See  Appendix  A) 


(6) 


where  a2,  A2&  represent  the  hies  errors  in  verticil  velocity  and  acceleration. 
The  first  two  terms  of  Equation  6  are  the  noise  error  terms  and  the  third 
term,  the  bias  error  term. 

a)  Disc  ass  ion  oi  Variables  in  Equation  6 

The  noise  error  in  vertical  velocity  (r*  )  and  acceleration  (<r~) 
depend  upon  the  noise  present  in  die  radar  coordinates  (e  ),  the  type  of  smooth¬ 
ing  technique  used  (polynomial  of  degree  P),  the  number  of  data  points  used  in 
the  smoothing  process  (N),  and  the  time  spacing  between  consecutive  data  points 
(At).  Of  these  four  variables,  two  are  dependent  upon  the  radar  being  used. 

For  example,  for  an  FPS  16  radar,  v  varies  between  10  and  15  meters,  depend¬ 
ing  upon  slant  range  and  tracking  conSitions,  and  At  is  generally  fixed  at  1/10 
of  a  second.  The  other  two  variables,  the  degree  polynomial  and  the  number 
of  data  points  are  arbitrary  and  will  be  chosen  in  such  a  way  so  as  to  minimise 
Equation  6.  The  values  of  Z  and  Z  in  Equation  6  depend  primarily  upon  the 
sphere  apogee,  but  are  also  affected  to  a  much  lesser  degree  by  atmospheric 
density.  If  a  balloon  flight  in.the  1962  Standard  Atmosphere  from  a  given 
apogee  is  assumed,  then  the  Z  and  Z  computed  from  this  theoretical  flight 
will  be  very  nearly  the  same  as  that  of  an  operational  flight  having  the  same 
apogee,  since  deviations  of  the  ac.  '.tmospheric  profile  from  the  1962 
standard  will  have  a  negligible  effect  on  sphere  velocity  and  acceleration,  (the 
notation  will  refer  to  the  density  profile  of  the  1962  Standard  Atmosphere). 


The  bias  errors  in  velocity  and  acceleration  (a£,  u  t)  depend  upon 
N.  At,  the  degree  polynomial,  and  the  true  position  Z  vs  t  of  the 
sphere.  The  true  position  path  is  primarily  a  function  of  sphere  apogee,  and 
only  secondarily  dependent  upon  the  deviations  of  density  from  a  standard  atmos. 
Since  the  true  position  field  of  a  passive  sphere  is  never  known,  some 
assumption  must  be  made  so  that  the  bias  error  term  can  be  evaluated.  A 
representative  position  field  can  be  obtained  for  the  sphere  assuming  that  the 
sphere  fell  in  the  1962  Standard  Atmosphere  from  a  given  apogee.  Using  this 
assumption,  the  next  section  will  describe  the  technique  for  evaluating  the 
bias  error  in  density. 

b)  Estimation  of  Noise  and  Bias  Error  Terms 

Equation  6  has  been  determined  to  be  a  function  of  N,  At,  degree 
polynomial,  and  apogee.  The  problem  to  be  solved  is:  determine  the  N.degree, 
and  apogee  that  will  minimize  Equation  6  in  the  70-100  km  region,  assuming  a 
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known  radar  error  and  a  known  sample  rate  (At).  For  the  passive  sphere 
system,  an  FFS  16  tracking  radar  is  assumed,  with  RMS  error  of  15  meters 
in  the  Z  component,  and  a  sample  rate  of  2  per  second  (At  =  1/2  second).  In 

t  actuality,  FPS  16  data  is  generally  digitised  in  1/10  second  intervals.  However, 

to  obtain  independent  measurements  and  to  minimise  the  data  handling  problem, 
five  1/10  second  data  points  are  averaged,  producing  independent  1/2  second 
data  points. 

Estimation  of  Noise  Error  Terms 

There  are  two  methods  of  estimating  the  noise  error  terms:  a)examining 
an  actual  flight  of  a  passive  sphere  tracked  by  two  identical  FPS  16  radars  and 
b)  using  equations  which  directly  relate  <r.  and  c-  to  N,  At,  <r  ,  and  the  smooth- 

Z  Z  2 

ing  function.  In  method  a9  densities  are  calculated  from  each  set  of  radar  data 
using  the  same  number  of  points  (N),  and  degree  of  smoothing  polynomial  (P). 
Since  the  same  bias  appears  in  the  density  computations  from  each  of  the  two 
radar  tracks,  calculation  of  the  RMS  difference  between  the  densities  as 
obtained  by  the  first  radar  and  those  from  the  second  radar  determines  the  noise 
error  terms.  The  equations  of  method  b  are  presented  in  Figure  3  and  derived 
in  Appendix  B.  As  seen  in  Figure  3,  the  form  of  the  noise  error  expression 
is  determined  by  the  degree  (P)  of  the  smoothing  polynomials  used.  Linear 
polynomial  smoothing  is  defined  as  fitting  a  linear  polynomial  over  N  data  points 
and  assigning  the  slope  of  the  fit  to  be  the  velocity  at  the  midpoint ^N+l  ^  of  the 

«  interval.  Linear -linear  smoothing  to  obtain  acceleration  is  described  as 

fitting  N  position  points  to  a  linear  polynomial  to  obtain  velocities  and  fitting  M 
of  these  velocities  by  a  linear  polynomial  to  obtain  acceleration  (see  Ref.  1). 

A  cubic -linear  fit  is  described  as  fitting  N  position  points  to  a  cubic  polynomial, 
taking  the  slope  at  the  midpoint  as  the  velocity,  and  fitting  M  of  these  velocities 
by  a  linear  polynomial  to  obtain  acceleration.  The  method  is  similar  for  cubic - 
cubic  and  linear -cubic  smoothing  techniques.  (Second  and  fourth  degree  poly¬ 
nomial  smoothing  result  in  the  same  slope  at  the  midpoint  of  the  interval  as 
does  the  first  and  third  degree  smoothing,  respectively. )  Quadratic  smoothing 
is  defined  as  fitting  a  second  degree  polynomial  to  position  points,  ar.d 
evaluating  the  first  and  second  derivatives  of  the  polynomial  (at  the  midpoint) 
as  the  velocity  and  acceleration,  respectively.  The  validity  of  these  formulas 
has  been  established  by  comparisons  with  RMS  errors  which  were  obtained  by 
method  a. 

Estimation  of  Bias  Errors 

For  a  given  apogee,  bias  errors  can  be  determined  by  the  following 
technique  (see  Figure  4).  Given  a  drag  table  and  sphere  apogee,  the  equations 
of  motion  can  be  integrated  to  obtain  the  theoretical  path  of  the  sphere,  assuming 
1962  Standard  Atmosphere.  The  vertical  position  Z,  velocity  Z,  and  acceleration 
2  are  determined  by  the  theoretical  trajectory  program.  The  Z  position 
coordinates  (Z,t)  are  now  tre  \ted  as  radar  data  and  the  smoothing  routine  is 
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Figure  3:  Noise  Error  Formulas  for  Polynomial  Smoothing.  (Derivations 
Given  in  Appendix  B) 
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applied  using  N  data  points  and  the  degree  polynomial  F.  The  smoothed  Z,  Z, 
and  Z  coordinates  differ  from  the  theoretical  Z,  Z,and  Z  coordinates  only 
because  of  the  bias  error  resulting  from  the  smoothing  technique.  (No  noise 
has  been  introduced  into  the  data. )  The  smoothed  coordinates  are  then  re -sub* 
stituted  into  the  equation  of  motion,  using  the  same  drag  table,  and  density  is 
computed.  The  only  difference  between  this  computed  density  and  the  original 
input  density  (1962  Standard  Atmosphere  Density)  is  due  to  the  bias  error 
induced  by  the  smoothing  function.  The  percent  bias  error  in  density  is  then 
determined  as  the  ratio  of  the  computed  minus  the  standard  density  to  the 
standard  density. 


c)  Determination  of  Optimum  Smoothing  Technique 

The  optimum  smoothing  for  determining  density  is  defined  to  be 
that  smoothing  of  vertical  coordinates  that  minimizes  the  sum  of  the  noise  and 
bias  error  variances  in  density  (i.  e.  minimizes  Equation  6).  The  determination 
of  the  polynomial  and  interval  that  minimizes  Equation  6  in  the  70-100  km  region 
was  derived  by  employing  the  above  techniques.  Initially,  an  escape  altitude 
of  125  km  was  chosen.  For  each  type  double  smoothing  (linear -linear,  cubic  - 
linear,  linear -cubic,  cubic -cubic),  the  total  error  in  density  due  to  velocity  and 
acceleration  errors  (Equation  6)  was  computed  for  all  possible  combinations  of 
N  and  M.  The  noise  error  was  calculated  by  method  b  (the  formulas),  as  given 
in  Figure  3,  using  At  =  1/2  second  and  <r  =15  meters.  The  bias  error  was 
computed  as  described  by  Figure  4.  Figures  5-8  are  examples  of  the  total 
error  plots  that  were  generated.  The  plots  indicate  the  percent  bias  error  (the 
deviation  of  density  ratio  from  1)  and  the  l<r  confidence  bands  of  the  noise  error 
about  the  bias.  Careful  analysis  of  total  error  plots  for  all  combinations  of 
degree  polynomial  (double  smoothing)  and  N  and  M  resulted  in  the  choice  of  the 
19-21  linear -cubic  combination  as  optimum. 

Total  error  plots  were  then  generated  in  the  same  fashion,  using  a 
quadratic  polynomial  fit  and  its  first  and  second  derivatives  for  velocity  and 
acceleration.  The  best  smoothing  interval  for  using  a  quadratic  was  determined 
to  be  31  data  points  (Figure  9).  In  comparing  the  optimum  quadratic  and  the 
optimum  linear -cubic  smoothing  techniques,  it  is  easily  seen  that  the  19-21  linear - 
cubic  produced  significantly  better  results  in  the  70-100  km  region.  The  probable 
explanation  for  this  is  the  following.  By  fitting  two  different  functions,  one  to 
get  velocity  and  the  other  to  get  acceleration,  it  is  possible  to  partially  compen¬ 
sate  for,  say,  a  positive  bias  in  density  due  to  a  velocity  error,  by  using  a 
different  degree  polynomial  or  different  interval  to  generate  accelerations  which 
will  produce  a  negative  density  bias.  This  advantage  is  not  present  when  using 
a  single  function  for  smoothing. 

The  decreased  density  accuracy  below  60  km  for  the  19-21  linear- 
cubic  is  due  to  the  slower  vertical  velocity  of  the  sphere  which  results  in  the 
constant  time  smoothing  technique  applied  over  a  shallower  layer  of  altitude. 
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Figure  5:  Blue  end  Noise  Errors  in  Density  for  19-21  Linear -Cubic  Smoothing 


Figure  7:  Bias  and  Noise  Errors  in  Density  for  31-7  Linear-Linear  Smoothing 
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Figure  8:  Biu  and  Noise  Errors  in  Density  for  49-11  Cubic-Cubic  Smoothing 
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Figure  9:  Bias  and  Noise  Errors  in  Density  for  31  Point  Quadratic  Smoothing 


This  situation  was  alleviated  by  expanding  the  position  smoothing  interval  from 
19  points  at  all  altitudes  above  60  km  to  a  maximum  of  51  points  at  approximately 
35  km  in  such  a  fashion  as  to  maintain  a  2%  noise  error  in  density. 

The  19-21  linear -cubic  smoothing  interval  was  next  evaluated  for  other 
sphere  apogees  from  70  to  150  km,  and  the  noise  and  bias  plots  were  generated. 
As  anticipated,  an  increase  in  sphere  apogee  resulted  in  higher  velocities  and 
accelerations  experienced  by  the  sphere  at  a  given  altitude,  which  increased 
the  bias  error  in  a  measured  density.  Conversely,  from  observing  the  noise 
error  terms,  it  is  seen  that  an  increase  in  velocity  and  acceleration  results 
in  a  decrease  of  the  noise  error.  In  calculating  the  total  error  in  the  70  to  100 
km  region  for  the  various  apogees,  it  was  determined  that  the  125  km  apogee  is 
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most  practical  (See  Table  4).  Not  only  is  the  total  error  one  of  the  smallest 
at  every  altitude,  but  from  the  logistic  standpoint,  the  125  km  apogee  -S  more 
easily  obtained  by  available  rocket  hardware  systems  than  higher  apogees. 
The  actual  noise  and  Was  errors  in  density  for  the  various  apogees  are  found 
in  Table  4.  This  data  verifies  the  requirement  that  serious  degradation  of 
density  accuracy  does  not  result  for  other  apogees. 


Tabl*  4 

Percent  Noisa,  BU»,  end  ToUl  Error  in  Density  for  Sphere  Apogee*  from  70  to  ISO  Km. 


■ 

.400.  i 

: sr~ 

...  ?i __ 

Total 

Sift* 

90  _ 

r: 

»s_  _ 

'  Total 

Bias 

.  _  _.«0_  1 

T0UI 

75 

ToUl 

BUS 

70 

Nu°fM  Tefal 
Bla* 

_  <10 _ 

NoiM  Total 

Bias 

ISO 

— 
m  m 

♦Is  4 

.1  • 

«.J  «•> 

Vs  > 

ail 

♦!  15 

QB 

140 

on 

m 

SB 

Dl 

♦U  2 

DDI 

12S 

HD 

,u  • 

DO 

.5  * 

m 

Em 

l 

m 

i  • 

US 

'A  » 

HD 

as 

♦I  s 

HD 

QB 

EH 

too 

—  - 

:: 

.?  * 

»  2.S 

42 

a  t 

HD 

IB 

i  * 

90 

mm 

HD 

»•*  1.  s 

>•*  1  5 

♦  .5  *'5 

m 

SO 

HD 

is  .  . 

♦  .S  *‘5 

!  1 

70 

:  » 

f 

! 

•* 

! 


Total  Error 


*e*  +  Si***] 


/* 


2.  2.  4  Vertical  Wind  Error  Term 


The  density  error  variance  resulting  from  vertical  winds  is  given 
by  the  expression 


(See  Appendix  A) 


(7) 


To  a  falling  sphere,  a  vertical  wind  appears  identical  to  a  change  in 
density.  As  a  result,  a  data  reduction  program  cannot  distinguish  density 
perturbations  from  vertical  wind  oscillations.  In  order  to  compute  densities, 
an  assumption  must  be  made  concerning  either  vertical  winds  or  density  per¬ 
turbations  in  the  atmosphere. 
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Assumption  A 

Assuming  no  vertical  motions  in  the  atmosphere,  Equation  1  can  be  solved 
by  substituting  W  =0  on  the  right  side  of  Equation  l  and  evaluating  all  other  terms 
by  conventional  means.  Under  this  assumption,  any  vertical  winds  present  in 
the  atmosphere  will  appear  as  density  perturbations.  The  relationship  between 
vertical  winds  (amplitude  and  wavelength)  and  density  perturbations  is  exhibited 
in  Figure  10  for  a  sphere  apogee  of  125  km.  Figure  10  also  contains  the  effect 
of  the  19 -21  (expanding  to  51-21)  linear -cubic  smoothing  filter  on  the  vertical 
wind.  An  example  will  clarify  this  effect.  If  a  sinusoidal  vertical  wind  of  wave¬ 
length  2  km  with  amplitude  of  1  m/sec  is  present  at  50  km,  then  this  vertical 
wind  would  be  damped  by  the  smoothing  and  appear  as  having  amplitude  0.  52 
m/sec.  (obtained  from  Figure  11).  Since  the  program  attributed  the  0. 52  m/ sec. 
vertical  motion  as  a  density  perturbation,  the  result  of  the  actual  l  m/sec. 
vertical  wind  would  be,  using  Equation  7,  a  l.  3%  density  perturbation  as  shown 
in  Figure  10.  The  density  error  introduced  by  a  vertical  wind  of  amplitude 
X  m/sec.  is  obtained  by  multiplying  the  error  introduced  by  a  l  mf  sec.  vertical 
wind  by  X. 


vVavc  length  km 


Figure  10:  Density  Error  Produced  by  a  Vertical  Wind  of  1  m/sec  -  125  km  Apogee 


Assumption  B 

If  density  is  assumed  to  follow  some  mean  path,  then  perturbations 
from  this  path  can  be  attributed  to  vertical  winds.  Since  density  varies 
exponentially  with  altitude,  a  mean  exponential  path  is  appropriate. 
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Using  this  assumption,  vertical  winds  can  be  computed  by1  the  equation 


Z  -  ± 

w  -  ma*a 

z  •  cl  z 
?KZpo* 

A  description  of  the  variables  in  this  equation  and  its  application  is  given  in 
Reference  1. 


Since,  according  to  present  knowledge,  meteorologists  accept  density 
perturbations  at  least  as  much  as  they  accept  vertical  winds,  Assumption  A, 
that  is  Wz  =  0y  has  been  incorporated  into  this  program.  The  magnitude  of 
vertical  winds  in  the  upper  atmosphere  is  not  yet  well  substantiated.  Consequently, 
significant  density  errors  may  result  when  using  this  assumption  if  the  vertical 
winds  exceed  a  few  m/sec  (See  Figure  10). 

2.  3  Summary  of  Density  Errors  Using  Optimum  Smoothing 

The  total  percent  error  in  density  resulting  from  a  computation  using 
the  High  Altitude  ROBIN  program  with  optimum  smoothing  cannot  be  precisely 
determined  because  of  the  unknown  accuracy  of  the  drag  table  and  the  occurrence 
of  unknown  vertical  winds.  The  drag  table  has  been  taken  from  the  work  of 
Heinrich  in  the  supersonic  mach  number  regime,  and  from  Lawrence  and  Goin 
in  the  subsonic,  with  values  interpolated  and  extrapolated  between  these  two  drag 
tables  when  necessary.  The  accuracies  specified  by  the  experimenters  are: 
Subsonic,  ±2%  RMS  error;  and  supersonic  maximum  possible  error,  from  ±2.  3% 
to  ±27.  9%,  depending  on  Reynolds  number. 


The  vertical  wind  error  depends  upon  the  magnitude  of  the  vertical  wind 
and  the  altitude  at  which  it  appears.  A  2  m/sec. ,  10  km  vertical  wind  will 
result  in  a  density  error  less  than  2%  above  70  km;  however,  at  an  altitude  of 
40  km,  the  same  2  m/sec.  vertical  wind  will  result  in  a  9%  error  in  density. 

The  error  in  density  due  to  errors  in  horizontal  winds  and  velocity  is  less 
than  1%  at  all  altitudes. 

The  error  in  density  due  to  a  10%  error  in  the  initial  estimate  of  temperature 
becomes  less  than  1%  after  4  km  of  flight  data. 

The  other  contributing  terms  to  density  errors  have  been  accurately 
determined  for  the  19-21  linear -cubic  (which  expands  to  a  maximum  51-21  at 
35  km)  smoothing.  The  noise  error  in  density  is  less  than  2%  to  90  km  and 
less  than  8%  to  100  km,  and  the  bias  error  is  less  than  1%  to  70  km  and  a  vary¬ 
ing  bias  from  lto5%  from  70  to  100  km.  The  amount  of  detail  in  the  density 


that  passes  through  the  linear«cubic  filter  is  determined  by  the  frequency 
response  curves  (Figure  11).  For  a  given  wavelength  of  a  density  oscillation. 
Figure  11  can  be  used  to  determine  the  percent  of  the  amplitude  of  the  oscillation 
that  will  filter  through  the  reduced  data  as  a  function  of  altitude.  As  seen 
i  n  Figure  11,  wavelengths  less  than  10  km  are  largely  damped  by  the  smooth* 
ing  at  altitudes  exceeding  70  km,  whereas  below  70  km  wavelengths  of  the  order 
of  S  km  are  easily  observed. 
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2.  4  Density  Validity:  Sphere  Collapse  Checks 

Density  data  is  defined  to  be  valid  if  the  sensing  sphere  is  properly 
inflated.  A  collapseelor  elongated  sphere  will  result  in  highly  ficticious  density 
measurements.  Two  checks  are  used  to  determine  proper  sphere  inflation: 

A  time  of  fall  check  and  a  density  gradient  check  (Xcheck). 

2.  4. 1  Time  of  Fall  Check 
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Experience  with  the  Areas  ROBIN  System  has  shown  that  a  collapsed 
sphere  falls  with  a  noticeably  slov'er  vertical  velocity  than  does  an  inflated  sphere. 
As  a  result,  it  is  generally  possible  to  determine  sphere  collapse  by  comparing 
the  fall  velocity  of  an  actual  sphere  to  the  estimated  fall  velocity  for  an  inflated 
sphere.  Estimated  fall  velocities  for  the  ROBIN  sphere  have  been  obtained 
from  computer  simulated  flights  assuming  the  1962  Standard  Atmosphere  and 
±20%  deviations  in  density  from  the  1962  Standard  Atmosphere.  The  times  of 
fall  for  sphere  apogees  between  150  and  70  km  is  given  in  Table  5.  Above  80  km 
the  time  of  fall  is  strongly  dependent  upon  the  sphere  apogee.  In  certain  cases 
an  error  of  a  few  km  in  determining  sphere  apogee  could  change  the  estimated 
time  of  fall  more  than*20%  deviation  in  density  from  the  1962  Standard  Atmosphere. 
For  this  reason,  the  time  of  faU  test  is  not  used  above  80  km.  The  altitudes 
where  the  time  of  fall  test  is  used  for  each  sphere  apogee  is  indicated  by  the 
heavy  lines  of  Table  5.  If  the  time  of  fall  of  an  actual  sphere  does  not  lie 
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between  thxt  estimated  for  a  sphere  falling  through  an  atmosphere  whose  density 
is  0. 8  Ph£2"  an<*  *  density  of  1. 2  p,,^  „  sphere  collapse  is  indicated. 

The  performance  of  the  time  of  fall  test  has  been  evaluated  for  numerous 
flights  with  varied  apogees.  In  Figure  12  a  comparison  is  made  between  collapse 
as  indicated  by  the  time  of  fall  test,  the  X  check  ( discussed  in  section  2. 4.  2), 
and  observations  of  the  density  ratio.  Experience  has  shown  that  a  sharp  increase 
in  the  density  ratio  over  a  few  kilometers  is  associated  with  sphere  collapse.  Of 
the  14  flights  that  were  investigated,  the  time  of  fall  test  accurately  determined 
collapse  on  12  flights.  On  the  other  two  flights  the  time  of  fall  test  was  apparently 
5  to  10  km  premature.  The  time  of  fall  test  can  only  pinpoint  collapse  to  within 
S  km.  If  collapse  occurs  near  the  bottom  of  the  test  layer  it  may  not  vary  the  time 
of  fall  sufficiently  to  detect  collapse  until  the  next  5  km  test  layer.  For  this 
reason  the  time  tolerances  which  indicate  collapse  are  quite  restrictive.  For 
spheres  that  never  inflated  or  spheres  whose  collapses  is  obvious  (by  the  three 
tests),  the  time  of  fall  greatly  exceeds  the  allowed  tolerances  for  inflated  spheres. 


TABLE  5 

Time  of  Fall  In  Seconds  Through  Standard  Atmoapbare 
and  220%  Deviation  from  Standard  for  Sphere  Apogees  Between  7C  and  150  Km 
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•  Indicates  collapse  by  time  of  fall  test. 
Holloman  Flights  are  ARCAS  ROBINS 


From  Figure  12  this  is  apparent  at  all  altitudes  and  is  a  strong  arguement  to 
support  the  validity  of  the  sphere  collapse  check.  A  debatable  point  is  whether 
the  tolerances  should  be  increased  slightly  at  the  expense  of  occasionally 
determining  collapse  5  km  too  late. 

2. 4. 2  Density  Gradient  Check  (X  check) 

Engler  (Ref.  1)  has  shown  that  below  60  km  the  density  gradient  can 
be  estimated  by  X  defined  as: 

x  ■  ttz -  ( ^  <*2>  -  ^  £■  <8> 

2  1 

An  estimate  of  the  value  of  X  one  would  anticipate  for  an  inflated 
sphere  is  determined  by  observing  the  density  gradient _1  dp  of  the  1962 

p  ds 

Standard  Atmosphere.  If  the  sphere  is  properly  inflated,  then  the  calculated  X 
should  approximate  this  density  gradient.  A  collapsed  balloon  will  result  in 
incorrect  density  data  and  consequently  X  will  deviate  considerably  frr  n  its 
anticipated  value.  To  allow  for  real  changes  in  the  density  gradient  ru  he 
atmosphere  from  the  1962  Standard,  the  change  in2_  dp  due  to  a  15%  perturb- 

p  dz 

ation  in  the  density  ratio  .  p  .  over  a  5  km  layer  (Figure  13)  has  been 

P„62» 

calculated  and  is  presented  in  Table  6. 


Figure  13:  Density  Profile  of  15%  Decrease  in  Density  Ratio  per  5  km 
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Table  6  shows  the  density  gradient  of  the  1962  Standard  and  the  change  occuring 
in  the  gradient  due  to  a  *  15%  perturbation  per  5  km. 


TABLE  6 

Density  Gradient  of  "62"  Standard  Atmosphere 
and  Changes  Produced  by  *15%  Density  Perturbation  Per  5  KM 


(}  *)  62  Std. 

gradient  for  -15%  gradient  for  +15% 

Perturbation  in  density  pnrturbntion  in  density 

difference 

60-55 

.00012 

.00009 

.06015 

6.00003 

55-50 

.00013 

.00010 

.00016 

6.00003 

SO -45 

.00014 

.00011 

.00017 

6.00003 

45-40 

.00015 

.00012 

.00018 

6.00003 

40-35 

.00916 

.00013 

.00019 

6.00003 

35-30 

.00015 

.00612 

.00018 

6.00003 

From  observation  of  Table  6  au  anticipated  value  of  X  =  .  00014  was 
chosen  and  a  tolerance  of  4J00005  allowed  for  density  perturbations  in  the  real 
atmosphere. 

In  addition  to  allowing  a  tolerance  in  X  for  real  density  perturbations 
a  tolerance  must  also  be  allowed  for  the  inaccuracies  of  the  parameters  used 
to  compute  X. 

It  can  be  shown  that  errors  in  a  calculation  of  Xare  primarily  a  result 
of  errors  in  the  velocity  measurements  (  Z^  and  2^)  rather  than  the  position 

measurements  ( and  Zj).* 


Considering  errors  in  X  as  a  function  of  velocity  errors,  the  error 
variance  formula  for  X  is  derived  as  follows 


ai  a-  x  a  • 

AX  =  -y—  Az,  +  — r~  Az, 

&Z2  2  3Zt  l 


Az. 


AX  = 


VZ1 


A\ 

t. 


(9) 


Since  the  program  is  designed  to  maintain  a  constant  percent  density 
error  { ££-)  of  .  02  below  60  km,  it  is  desirable  to  transform  Equation  9  into 
terms  of  Ap  _ 

P 


*  Zj  and  are  the  respective  first  and  last  points  of  the  position  smoothing  interval. 
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Assuming  p  =  ct  an  approximate  solution  to  the  equations  of 
motion  yields 

Z  =  (JL)1/2  e  p  1/2  (See  Ref.  1) 

kc  kc 

where  g,  k,  and  c  are  constants  defined  in  Reference  1. 

The  error  increment  in  £  is  approximated  as 

a-  3*  *  g  .1/2  Ap 

az  =  —  ap  =  (—j)  rth 


kc  2p" 


Thus 


-=  1/2 


Ao.. 


Substituting  (10)  into  Equation  9  yields 


AX,  =— — 
Z2-Zi 


Taking  the  variance  of  both  sides  of  Equation  11  and  recalling  that  Ap  is 
maintained  constant  below  60  km,  gives  ^ 


o-  2  =  — 1 _ 2  2(-2£-Y 

X  Tpi 


Thus 


where 


Az  =  Z_  -  Z  . 
2  1 


For  a-—- =.  02  the  noise  error  in  X  becomes 


02828 

Az 
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The  3 «•  noise  error  tolerance  for  X  ia 


.08484 
A z 


The  X check  ia  summarized  as  followa: 

Expected  value  of  X=  .  00014 

Tolerance  for  variations  in  density  ±  .  00005 

3<r  noise  tolerance  . .  08484 

. . — —  • 

Az 

A  sphere  is  considered  inflated  if 


.  00014  -  -/(.  0005)' 


. 08484 
Az 


<X<.  00014  +  /{.  00005 r 


,08484’ 

Az 


The  performance  of  the  X  check  is  evaluated  in  Figure  12.  For  9  of  the 
14  flights  the  collapse  altitude  as  determined  by  the  X  check  is  in  excellent 
agreement  with  the  collapse  altitude  as  indicated  by  time  of  fall  and  the 
density  ratio.  For  these  cases  the  X  check  appears  to  determine  the  collapse 
altitude  to  within  1  km  accuracy.  For  the  other  flights  the  X  check  generally 
indicates  collapse  several  kilometers  before  the  time  of  fall  check  .  This  is 
not  caused  by  radar  noise  but  is  undoubtedly  the  result  of  an  interesting 
meteorological  observation;  either  a  vertical  wind  or  density  perturbation. 

For  this  reason  the  X  check  should  be  interepreted  in  conjunction  with  results 
of  time  of  fall  test.  The  time  of  fall  test  is  the  primary  determination  of 
collapse.  The  X  check  serves  to  verify  the  performance  of  the  time  of  fall 
test  as  well  as  pointing  out  interesting  meteorological  phenomena. 
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Tbe  error  present  in  a  wind  computation  is  a  result  of  the  inaccuracies 
of  parameters  in  Equation  13.  The  inaccuracies  of  several  of  these  parameters 
have  a  negligible  effect  upon  the  wind  accuracy. 

3. 1  Variables  Not  Producing  Significant  Wind  Errors 

The  errors  in  computing  g  ,  the  gravitational  acceleration  ii^  the 
direction  perpendicular  to  the  tangent  plane,  is  less  than  .  0004  m/sec  .  This 
error  is  negligible  when  compared  to  the  larger  error  in  measuring  2 .  The 
error  in  g  is  approximately  1/10  the  magnitude  of  the  error  in  g  ,  and  for  the 
same  reason  can  be  ignored. 

The  Z  component  buoyancy  term,  p  V  g  /m  is  very  small  (approximately 
2  1 

.  006  m/sec  ).  The  error  in  computing  this  term  is  many  times  smaller  than  the 
term  itself  and  is  negligible  when  compared  to  other  larger  error  terms.  The 
error  in  the  X -component  buoyancy  term  p  g^/m,  is  approximately  1/10  the 

magnitude  of  the  error  in  the  term  p  g  / m,  and  similarly  can  be  ignored 
when  investigating  wind  accuracy.  J 

For  a  ROBIN  sphere  with  apogee  125  km,  the  coriolia  acceleration 
component  C  has  approximate  magnitude.  09  m/ sec^  at  100  km,  .  11  at  95  km, 
and  .  11  at  90  ^m.  Depending  upon  the  magnitude  of  the  denominator  of  Equation 
13,  C  can  make  a  large  contribution  to  wind  measurement.  If  ignored,  the 
coriofis  acceleration  at  100  km  would  result  in  a  ficticious  horizontal  wind  of 
magnitude  40  m/sec;  at  95  km,  a  ficticious  horizontal  wind  of  19  meters  per 
second;  and  at  90  km,  a  wind  of  10  m/sec.  Hence,  the  accuracy  of  C  must  be 
investigated.  is  computed  by  equation 

C  =  2  wlZcos  6  +  Y  sin  e)  (14) 

x  '  ' 

where  ^ 

W  =  rotation  rate  of  earth  =  7.  29  x  10  radians  per  sec. 

0  =  latitude 

The  error  in  C  is  a  result  of  the  error  in  Z.  For  example,  if  Z  were  in  error 
x 

by  10  m/sec,  the  resulting  error  in  C  would  be  only  .  0014.  Since  this  is  at 

•  • 

least  an  order  of  magnitude  less  than  the  error  in  X,  it  is  not  a  significant 

source  of  error  in  wind  calculation.  The  error  in  the  vertical  coriolis 

component  C  is  approximately  1/10  the  error  in  C  .  Consequently,  errors  in 
z  x 

both  coriolis  components  of  acceleration  are  not  significant  sources  of  wind 

measurement  error. 

The  effect  of  vertical  winds  on  horizontal  wind  measurements  depends 
upon  the  altitude  of  the  sphere.  Table  7  shows  the  magnitude  of  errors  that  will 
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result  from  a  5  m/eec  and  1  m/sec  vertical  wind  as  a  function  of  altitude. 

The  error*  may  get  as  large  as  a  couple  of  m/sec  at  100  km.  The  remaining 
variables  in  Equation  13  which  may  contribute  significantly  to  wind  errors  are 
the  horizontal  and  vertical  components  of  velocity  and  acceleration. 


tasle  ? 

XcprtHnlilin  Value*  ol  Horisental  Wind  Error 
Suuitiof  from  1  and  5  m/a* *  Vertical  Wind* 

»  5  m/**c 


Altitude 

km 

t 

i-Vk 

Representative  Vain* 
of  O  *ee  below 

Error  la  Vt 

too 

.500 

5 

.01 

250 

2. 5  m/eec 

90 

•625 

5 

.006 

200 

1.6 

to 

.550 

5 

.01 

160 

1.6 

10 

.250 

5 

.02 

100 

2.0 

60 

.160 

5 

.03 

5 

.15 

50 

-100 

5 

.05 

1 

.05 

40 

.  41 

5 

.12 

.05 

.007 

JO 

.  14 

5 

36 

.04 

.014 

w  * 
s 

1  m/eec 

Altitude 

km 

t 

W1 

Representative  Value 
of  O  see  below 

Error  in  W. 

100 

.500 

1 

.002 

250 

.5 

90 

•625 

1 

.0016 

200 

.32 

•0 

-550 

| 

.0011 

160 

.266 

70 

.250 

1 

.004 

100 

.4 

40 

.160 

1 

.0062 

5 

.03 

50 

.100 

1 

.01 

1 

.01 

40 

.  41 

1 

.  024 

.06 

.002 

JO 

.  14 

1 

.071 

.04 

.003 

EaLh. 


O  refers  to  the  term  m  ).  |  .  *■  ^  ^  is  the  percent  error 


-*  -  8s 

in  a  resulting  by  assuming  ■  0. 
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3.2 


Wind  Error  Equation 


« 


9 


Having  retained  only  those  error  components  which  are  the  predominant 
source  of  error  in  a  computed  wind,  the  wind  equation  simplifies  to 


W 

x 


=  X 


X  ± 

z-g 


The  first  order  of  approximation  to  an  error  in  W  resulting  from  the 
errors  in  the  parameters  is  given  by  Equation  15 


dw  s  dx  - 
x 


+ 


X  Z 

il 

(Z  -  g) 


2 


dz 


(15) 


If  the  differential  error  components  are  considered  as  random  noise  error  with 
normal  distribution,  then  by  taking  the  variance  of  Equation  15,  the  noise  error 
in  a  wind,  calculation  is  given  by  Equation  16  where  <r. ,  <r..  ,  etc.  are  the  noise 
errors  in  velocity  and  acceleration  due  to  the  noise  m  the  radar  data 


To  determine  the  bias  in  a  wind  measurement  resulting  from  the  biased  (over  - 
smoothing)  velocity  and  acceleration  measurements,  Equation  15  is  again  applied. 
Considering  the  component  differentials  as  the  bias  error,  the  square  of  the 
bias  wind  error  is  given  by  Equation  17. 


Az 


(17) 


where  the  Ax,  etc,  refer  to  the  bias  error  in  the  X  component.  The  total  wind 

error  (or  total)  is  defined  as  the  square  root  of  the  sum  of  the  noise  error 
w 

x 

variance  plus  the  bias  square  error  and  is  given  by  Equation  18 


<r  =  r 

W 

x  total  L 


<r 

w 

x 


2 


+ 


1/2 


(18) 
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The  optimization  problem  for  wind  determination  is  to  find  that  type  smoothing 
and  smoothing  intervals  for  horizontal  and  vertical  coordinates  that  minimize 
Equation  18.  As  in  the  case  of  the  density  smoothing,  the  noise  error  will 
decrease  as  the  smoothing  interval  increases  and  the  bias  error  increases  as  the 
smoothing  interval  increases,  so  that  a  minimum  does  exist  for  Equation  18. 

3.  3  Minimization  of  Total  Wind  Error 

3.  3. 1  Discussion  of  Variables  in  Noise  Error  Term 


The  noise  error  variance  (Equation  16)  is  a  function  of  the  vertical  velocity 
and  acceleration,  the  horizontal  velocity  and  acceleration,  and  the  errors  in 
measuring  the  horizontal  and  vertical  velocities  and  accelerations.  As  shown 
previously  (Figure  3),  the  noise  error  in  a  velocity  or  acceleration  measurement 
is  a  function  of  the  degree  of  smoothing,  smoothing  interval,  time  spacing  between 
data  points,  and  the  error  in  the  radar  coordinates.  The  types  of  smoothing 
considered  were  linear -linear,  cubic -linear,  cubic -cubic,  quadratic  and  quartic 
(quadratic  and  quartic  smoothing  consists  of  fitting  a  second  and  fourth  degree 
polynomial  respectively  to  the  position  coordinates  and  evaluating  the  first  and 
second  derivatives  of  the  fits  as  velocity  and  acceleration).  The  noise  error  formulas 
for  the  various  degrees  of  smoothing  for  velocity  and  acceleration  are  given  in 
Figure  3.  For  any  actual  sphere  flight,  and  for  any  of  the  smoothing  techniques 
considered  for  measuring  velocity  and  acceleration,  the  noise  error  in  winds 
(Equation  16)  can  be  calculated  by  knowing  the  radar  accuracy  fx  ,  x  *<r  )$time  spacing 
between  data  points  <1/2  second  for  the  independence  requirements ),^ancl  the 
velocity  and  acceleration  experienced  by  the  sphere  (these  can  be  taken  merely  as 
the  smoothed  velocity  and  acceleration  measurements).  Therefore,  by  taking  a 
particular  ROBIN  flight  and  knowing  the  radar  accuracies,  the  noise  error  variance 
can  be  calculated  for  any  length  velocity  and  acceleration  smoothing  intervals,  as 
well  as  for  linear,  cubic,  quadratic,  and  quartic  velocity  and  acceleration  smoothing 
procedures. 


3.  3.  Z  Discussion  of  Variables  in  Bias  Error  Term 


The  bias  error  term  has  been  given  by  Equation  17.  To  evaluate  this 
expression,  equations  must  be  derived  for  the  bias  error  in  the  velocity  and 
acceleration  measurements  (Ak,  Ax,  Az,  Az).  The  rest  of  the  variables  in  Equation 
17  are  immediately  available  from  the  smoothed  velocity  and  acceleration  values. 

The  bias  velocity  and  acceleration  measurements  can  be  written  as  a  function  of  the 
degree,  interval,  At,  and  the  true  position  field  of  the  sensor  versus  time.  If  one 
assumes  the  true  position  field  of  the  sensor  to  be  a  polynomial  of  degree  3,  and 
assuming  a  linear  polynomial  smoothing  to  evaluate  the  derivative,  then  Ax  defined 


*A  a  of  10  to  15  meters  in  each  component  is  a.  representative  estimate  of  FPS-16  radar 
accuracy  obtained  over  a  period  of  years  from  numerous  dual-tracked  flights  flown 
at  various  ranges. 


as  &  -  X  .  .  Is  given  by  A  At2  (3N2-7)/20  where  A,  is  the  coefficient 

of  thertl?ird  degree  term.  (See  Appendix  C).  If  cubic  smoothing  is  used,  then 
there  will  be  no  bias  error  for  position  fields  up  tu  and  including  fourth  degree 
polynomials.  That  is,  X  fmoothed  e^ua^  &  true  (a8#um*ng  ooiae  free  data). 
The  bias  error  formulas  for  single  and  double  fit  smoothing  are  presented  as 
Figure  14  and  derived  in  Appendix  C.  To  apply  the  bias  error  formulas  one  must 
be  able  to  substantiate  the  degree  position  field  and  determine  the  value  of  A^  , 
the  cubic  coefficient.  For  quadratic  and  quartic  smoothing,  the  coefficient 
A^  must  be  determined.  The  next  section  will  show  how  the  position  field  and  A^ 
can  be  estimated  for  linear -linear,  cubic -linear,  and  cubic -cubic  smoothing. 

A  similar  approach  has  been  employed  for  single  fit  smoothing  (quadratic  and 
quartic).  , 

3.  3.  3  Procedure  for  Evaluating  <r 

Wx  total 

The  total  wind  error  (Equation  18)  was  evaluated  for  linear -linear,  cubic - 
linear,  and  cubic -cubic  smoothing  with  N  (velocity  smoothing  interval)  and  M 
(acceleration  smoothing  interval)  varying  over  a  wide  range  of  values  to  determine 
the  optimum  combination  of  degrees  smoothing  and  smoothing  interval  that 
minimizes  Equation  18.  Linear -cubic  smoothing  was  not  considered  for  the 
following  reason.  Assuming  position  to  be  an  kth  degree  polynomial,  then  velocity 
is  an  k-lst  degree  polynomial  and  acceleration  k-2nd  degree  polynomial.  Since 
acceleration  is  of  one  degree  less  than  velocity,  one  does  not  anticipate  that  the 
degree  smoothing  on  acceleration  should  exceed  the  degree  smoothing  on  velocity, 
thus,  linear -cubic  is  not  considered. 

This  logic  is  not  applicable  for  the  case  of  smoothing  for  density.  For  that 
case,  linear -cubic  was  considered  because  the  problem  there  was  to  determine 
two  degrees  for  double  smoothing  that  would  effectively  cancel  the  bias  error  in 
density. 


Figure  15  shows  the  steps  in  logic  used  in  evaluating  the  total  error 
equation  for  winds.  The  discussion  will  be  restricted  to  the  X  coordinate  only, 
However  the  same  procedure  holds  for  the  Y  and  Z  coordinates.  Step  1  is  to 
determine  the  radar  noise  cr  ,  and  the  time  spacing  (At)  for  independent  data  points. 
These  two  quantities  will  permit  the  evaluation  of  the  noise  error  formulas  for 
linear  and  cubic  smoothing  and  for  any  N  and  M.  <r  is  most  easily  determined 
from  dual  tracked  flights  by  calculating  the  standard  deviation  of  the  differences 
in  the  X  position  values  from  the  two  tracks.  If,  however,  the  flight  being 
analyzed  is  tracked  by  only  one  FPS-16  radar,  then  theoretical  radar  position 
errors  of  perhaps  a  tenth  of  a  mil  in  azimuth  and  elevation  angles  and  4.  5  yards 
in  slant  range  may  be  assumed,  and  these  values  can  be  converted  to  rectangular 
coordinate  errors.  The  remaining  steps  of  Figure  15  basically  exhibit  the  logic 
and  procedures  required  to  evaluate  the  bias  error  components,  Ax  Ax, for  each 
type  double  smoothing.  In  order  to  evaluate  the  bias  errors  an  assumption  must 
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DOUBLE  SMOOTHING 


Position  assumption  for  N  <  N*  data  points 


X  =  A0  +  Ajt  +  A^t2  +  A^t3 


Velocity  Bias 
Linear  fit 


Ai  (3H2-7) 

1  20 


Ai3  =  0 


Cubic  fit 

Position  assumption  Quintic  for  M  data  points  where  2M  +  N  -  2  £N 

*  2  3  4 

X  =  Bq  +  Bjt  +  B^t  +  B^t  +  B^t  (velocity) 


Acceleration  Bias 


Linear -linear  fit 

**11**33 -*11 

Cubic -linear  fit 

AStj,  =  BjAt,2  <3M2-7> 
20 

Cubic -cubic  fit 

A*33  =  0 

SINGLE  SMOOTHING 

Position  assumption  for  N  data  points 

X 

*  A.  +  A.  +  A,t2  +  At3  +  A  .t4 

0  12  3  4 

Velocity  Bias 

Quadratic  fit 

Ai  -  A3At2(3N2.7) 

2  20 

Quartic  fit 

o 

ii 

4 

Acceleration  Bias 

Quadratic  fit 

A~  A  Aa2  (3N2-13) 

Ax,  =  A  .At  - - - 

2  4  7 

Quartic  fit 

© 

it 

4 

Figure  14:  Bias  Error  Formulas  for  Polynomial  Smoothing 
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Figure  15:  Schematic  for  Calculating  Total  Error  in  Winds 


be  made  concerning  the  path  of  the  sphere.  Step  two  is  to  assume  that  the  true 
position  path  of  the  sensor  is  a  third  degree  polynomial.  Step  three  is  to 
determine  the  maximum  number  of  points  (N)  for  which  this  assumption  is  valid. 

This  number  can  be  determined  by  comparing  the  theoretical  or  dual  track 
noise  error  («r  )  to  the  standard  deviation  of  the  position  data  from  the  cubic  fit 
over  N  data  points.  If  the  cubic  assumption  is  valid  for  N  data  points,  then  the 
noise  error  should  approximate  the  standard  deviation  of  the  position  points  from 
the  cubic  fit.  N  in  turn  is  increased  until  the  standard  deviation  of  the  position 
points  from  the  cubic  fit  significantly  exceeds  «r  .  The  smallest  value  of  N  for  which 
the  cubic  assumption  is  found  valid  over  all  segments  of  the  flight  is  called  N*. 

Step  four  is  to  choose  some  N  less  than  or  equal  to  N*  for  an  initial  velocity 
smoothing  interval.  Step  five  is  to  fit  a  cubic  polynomial  to  the  N  position  data 
points.  The  coefficients  of  the  cubic  fit  should  be  given  with  respect  to  t  =  o 
as  the  midpoint  of  the  smoothing  interval  in  order  to  utilise  the  bias  error 
formulas  of  Figure  19.  Thus,  the  velocity  at  the  midpoint  is  the  coefficient 
(step  6).  Since  the  cubic  polynomial  was  fit  to  assumed  cubic  position  points,  no 
bias  error  will  result  in  the  computed  velocity  A^.  The  error  in  the  cubic  velocity 
due  to  radar  noise  (<r  )  is  calculated  by  the  noise  error  formulas  in  Figure  3. 

W3 

The  noise  and  bias  errors  are  next  calculated  for  linear  velocity  smoothing  by 
steps  7  and  8.  Step  7:  A  linear  polynomial  is  fitted  to  the  cubic  position  data 
points.  In  this  case  there  is  both  a  noise  and  a  bias  error  associated  with  the 
linear  velocity.  Step  8:  The  noise  error  is  computed  from  the  linear  noise 
error  formula  (Figure  3).  The  bias  error  is  estimated  from  the  formula 

2  2 

Akj  =  At  (3N  -7)/20  (  Figure  14  )  where  A^  is  the  estimate  of  the  true  cubic 

coefficient  of  the  position  field  determined  by  fitting  a  cubic  polynomial  to  the 
assumed  cubic  position  data  which  is  contaminated  by  radar  noise.  Due  to  the 
noise  in  the  data,  the  A  used  to  calculate  the  bias  error  is  not  the  true  A 
of  the  position  field.  Therefore  there  will  be  an  error  in  the  calculation  of  the 
bias  error.  However,  the  variance  of  the  error  will  be  less  than  the  cubic  noise 
error  term  (Comment  Step  8).  Consequently,  even  though  the  bias  error  calculated 
is  only  an  estimate  of  the  true  bias  error,  bounds  can  be  placed  on  the  accuracy 
of  the  estimate.  Furthermore,  when  an  actual  flight  is  analyzed  and  a  large  number 
of  calculations  of  have  been  made  the  calculated  values  of  A^  will  have  as 
their  mean  the  true  A  so  that  the  calculated  bias  error  will  be  centered  about 
the  true  bias  error  ana  a  judgement  can  be  made  concerning  its  magnitude. 

For  the  sake  of  acceleration  error  calculations,  the  position  field  is  now 
assumed  to  be  of  degree  5.  Step  10  is  to  determine  the  maximum  number  of  points 
for  which  a  fifth  degree  polynomial  assumption  holds.  Call  this  number  N'.  N' 
is  found  by  the  same  technique  used  to  find  N*  (Refer  to  step  3).  M,  the  number 
of  data  points  used  to  compute  acceleration,  is  now  chosen  so  that  2M  +  N-2  is 
less  than  or  equal  to  N'  (Step  li).  The  number  2M  +  N-2  is  the  number  of  position 
points  that  influence  an  acceleration  calculation  when  M  points  are  used  for 
acceleration  calculations  and  N  points  for  velocity.  Position  assumed  quintic 
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over  N:  points  implies  velocity  is  quartic  over  N'  points.  Fitting  quartic 
velocity  by  a  cubic  polynomial  will  result  in  no  bias  error  (the  first  derivative 
of  a  cubic  polynomial  fit  coincides  with  the  first  derivative  of  a  quartic  poly¬ 
nomial  fit  when  evaluated  at  the  midpoint  of  the  interval ).  Step  12:  A  cubic 
polynomial  is  fitted  to  the  velocity  points  obtained  from  the  cubic  fits  of  position. 

No  bias  error  is  present  in  this  cubic -cubic  acceleration  (B^)  since  the 
velocity  points  fitted  were  free  of  bias  (Step  13 ).  The  noise  error  for  the  cubic  - 
cubic  acceleration  is  calculated  from  the  noise  error  formulas.  Next,  acceleration 
bias  using  cubic -linear  smoothing  is  determined.  Step  14:  A  linear  polynomial 
is  fitted  to  the  cubic  velocity  points  (i.  e. ,  velocities  obtained  by  cubic  fits  to 
position).  The  bias  error  for  the  cubic -linear  acceleration  is  given  by 

B  (3M^-7)/20t  where  is  the  estimate  of  the  true  cubic  coefficient  obtained 

by  fitting  a  cubic  polynomial  to  noisy  quartic  velocity  points  and  At^  is  the  time 
spacing  between  velocity  points,  (step  15).  Analagous  to  the  velocity  case  there 
will  be  an  error  in  the  calculation  of  the  bias  error,  since  the  calculated  is 
only  an  estimate  of  the  true  cubic  coefficient.  The  bounds  for  the  variance  of 
the  bias  estimate  can  be  shown  to  be  less  than  oA-  .  The  noise  error  for  cubic  - 
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linear  acceleration  is  calculated  from  Figure  3.  Step  16:  A  linear  polynomial 
is  fitted  to  the  velocities  obtained  by  the  linear  fits  of  position  data.  The  noise 
error  in  this  linear -linear  acceleration  is  given  by  Figure  3.  The  most  easily 
accessible  means  of  calculating  the  bias  error  for  the  linear -linear  acceleration 
is  to  difference  the  cubic -cubic  acceleration  from  the  linear -linear  acceleration 
(step  17).  Since  the  cubic -eubic  acceleration  contains  no  bias,  its  difference 
from  the  linear-linear  acceleration  will  estimate  the  bias.  The  error  variance 
bounds  of  the  estimate  of  the  bias  can  be  shown  to  be  less  than  a?.  .  Having 

X33 

calculated  the  noise  and  bias  errors  for  linear -linear,  cubic -cubic,  and  cubic - 
linear  velocity  and  acceleration  measurements,  a  substitution  for  these  values 
into  Equation  18  yields  the  total  wind  error  for  each  type  double  smoothing  (step  18). 
A  comparison  between  the  three  total  wind  errors  for  the  three  types  of  smoothing 
determines  that  type  smoothing  which  produces  the  minimum  total  error  for  the 
particular  N  and  M  chosen.  By  changing  the  values  of  N  and  M  the  procedure 
can  be  repeated  and  further  comparisons  made.  The  optimum  double  smoothing 
technique  is  that  combination  of  degree  polynomials  (cubic -cubic,  cubic -linear, 
linear -linear)  and  smoothing  intervals  (N  and  M)  that  gives  the  minimum  total 
wind  error.  Plots  of  the  total  wind  error  for  each  type  double  smoothing,  and 
for  N-M  values  of  51-35,  53-11,  and  31-11  are  presented  as  Figures  16-18.  These 
are  merely  throe  illustrations  of  some  of  the  possible  combinations.  Figures  19-20 
present  the  total  wind  error  for  quadratic  and  quartic  single  fit  smoothing.  The 
total  error  plots  for  single  fit  smoothing  are  determined  by  a  procedure  similar 
to  that  discussed  for  double  smoothing.  The  noise  and  bias  error  formulas  for 
single  fit  smoothing  are  also  given  in  Figures  3  and  14.  After  analyzing  plots 
of  the  types  illustrated  in  Figures  16-20  for  all  feasible  combinations  of  N  and  M, 
it  was  determined  that  the  51-35  cubic -cubic  provides  optimum  wind  reduction. 

With  this  type  smoothing,  the  total  wind  error  (which  consists  only  of  noise 
error)  remains  less  than  2  m/sec  to  85  km  and  10  m/sec  to  nearly  100  km.  Even 
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Figure  19:  Total  Error  in  X  Wind  and  Y  Wind  for  31  Point  Smoothing  of  Viper  Dart  13 
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Figure  20:  Total  Error  in  X  Wind  and  Y  Wind  for  99  Point  Smoothing  of  Viper  Dart  13 
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with  optimum  smoothing,  the  amount  of  detail  that  can  be  observed  at  the  very 
high  altitudes  is  limited  because  of  the  large  altitude  layer  used  in  the  smoothing. 
The  frequency  response  curves  (Figure  21)  indicate  the  detail  that  can  be  observed. 
Plotted  in  the  figure  is  the  ratio  of  the  amplitude  of  a  sinusoidal  wind  after  passing 
through  the  smoothing  filter  to  the  true  amplitude  of  the  original  wind.  Bach  curve 
giyes  the  ratio  versus  wavelength  at  a  specific  altitude.  For  example,  at  80  km 
the  amplitude  of  a  10  km  wave  in  the  wind  profile  would  appear  in  the  reduced  data 
to  have  i/4  its  actual  amplitude.  A  20  km  wave  would  retain  65%  of  its  amplitude 
in  the  reduced  data.  As  seen  from  the  figure,  at  altitudes  in  excess  of  70  kilo- 
meters,  wavelengths  less  than  10  km  are  largely  destroyed,  so  that  only  a  mean 
wind  profile  can  be  ascertained.  Below  70  km,  wavelengths  of  5  km  and  less 
will  appear  in  the  reduced  data. 

3. 4  Summary  of  Wind  Errors  Using  Optimum  Smoothing 

The  51-35  cubic -cubic  smoothing  technique  in  X,  Y,  and  Z  was  found  to  be  the 
optimum  for  measuring  wind  between  70  and  100  km.  Using  this  smoothing,  winds 
can  be  measured  to  at  least  an  accuracy  of  2  m / sec  to  85  km  and  10  m/sec  to  95  km. 
Above  70  km,  however,  fine  details  in  the  wind  profile  are  lost  so  that  only  details 
with  vertical  structure  in  excess  of  10  km  are  retained  by  the  smoothing  technique. 
Below  70  km,  fine  scale  structure  with  vertical  wavelengths  of  2  to  5  km  can 
be  observed  in  the  reduced  data. 
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Figure  21:  Steady  State  Wind  Error  and  Frequency 

Reeponee  for  51-35  Cubic-Cubic  Smoothing 
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Wind  Validity 


It  is  interesting  to  observe  that  a  collapsed  sphere  cannot  accurately 
measure  density  but  can  accurately  determine  wind.  The  reason  being  the 
following:  The  equation  used  to  compute  density  is  indirectly  proportional  to  the 
drag  coefficient.  When  a  sphere  collapses  its  drag  coefficient  changes  drastically. 
The  drag  coefficient  tables  available  are  for  spheres  only  and,  thus,  cannot  be 
used  for  density  determination  when  a  sphere  loses  its  spherical  shape.  On  the 
other  hand,  the  simplified  equation  (not  containing  buoyancy  and  coriolis  terms) 
used  to  compute  wind  does  not  depend  upon  the  drag  coefficient  of  the  sphere  and 
hence,  can  be  evaluated  even  when  the  sphere  is  in  a  collapsed  condition.  A  slight 
effect  of  sphere  collapse  on  wind  determination  does  occur,  however,  through 
the  buoyancy  term  which  requires  a  value  for  density.  The  insignificance  of  the 
buoyancy  force  above  30  km  discards  this  as  a  possible  source  of  wind  error. 

Thus,  wind  accuracy  can  be  maintained  for  a  sphere  in  collapsed  as  well  as 
inflated  condition. 


4.  PRESSURE 

The  equation  used  for  Pressure  Computation  is 

Z. 

pi  =  pi  i  +  r  1  pgd*  . 
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(19) 


In  order  to  use  Equation  19,  an  initial  pressure  (Po)  is  required.  This 
pressure  is  computed  from  the  gas  law 
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where  T  is  an  estimate  of  the  initial  temperature  generally  taken  from  the  1962 
Standard  Atmosphere  or  a  Supplemental  Atmosphere,  P  is  computed  from  the 
density  equation  (1)  by  using  T  to  evaluate  the  initial  Mach  and  Reynolds  numbers, 
R  is  the  universal  gas  constant  and  M  is  the  molecular  weight  of  air.  The 
integration  of  Equation  19is  accomplished  by  the  formula: 
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4.1  Pressure  Error  Equation 

In  order  to  simplify  the  derivation  of  the  error  equation  for  pressure, 
linear  integration  using  the  trapazoidal  rule  is  assumed,  i.  e. , 
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Thus 
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The  error  in  P  in  terms  of  density  errors  is 
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Two  sources  of  error  are  apparent  in  equation  20:  a)  the  error  in  pressure 
due  to  the  error  in  the  initial  computation  of  pressure  (from  the  initial  estimate  of 
temperature)  and  b)  the  error  in  pressure  due  to  the  integration  of  density  errors. 
Both  noise  and  bias  errors  are  present  in  density.  These  errors  will  effect  pressure 
differently. 

4. 1. 1  Pressure  Errors  Resulting  from  Error  in  Initial  Pressure 

As  seen  from  Equation  20,  the  error  in  pressure  at  any  point  (P.)  of  the 
flight  due  to  the  error  in  initial  pressure  is  constant.  The  percent  error  in  pressure 
due  to  the  initial  pressure  error  decreases  as  the  flight  progresses  since 


P. 
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becomes  very  small  after  a  few  kilometers.  For  a  10%  error  in  the  initial  estimate 
of  temperature  and  a  73  km  sphere  apogee.  Engler  (Ref.  1)  has  shown  that  after  6  to  7 
km  of  flight  the  resulting  error  in  pressure  has  diminished  to  less  than  1%.  For  a 
125  km  apogee,  theoretical  trajector-.es  have  been  computed  using  initial  temperatures 
of  Tq  =T„£2ii,  and  Tq  =  .9Tii^2»i»  and  T  =1.  lT^^at  100  km.  These  results  agree 
with  Engler's  and  are  presented  in  Table  8. 
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TABLE  • 


Percent  Error  In  Preeeure 
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4. 1.  2  Pressure  Errors  Resulting  from  Noise  Errors  in  Density 

Assuming  the  errors  in  each  density  data  point  (e  )  are  independent 

Pi 

with  mean  0  and  standard  deviation  o-  the  pressure  error  variance  at  a  point 

Pi 

P.  is  given  as 
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Equation  21  has  been  evaluated  for  the  expected  noise  error  profile  in 
density  for  the  expanding  19-21  linear -cubic  smoothing.  See  Table  4  for  pertinent 

It 

values  of  p  (100).  The  results  are  presented  in  terms  of  percent  pressure  error 
in  Table  9. 


TABLE  9 
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4. 1.  3  Pressure  Errors  Resulting  from  BicS  Errors  in  Density 

Bias  errors  arise  in  density  because  of  bias  introduced  by  the  smoothing 
function  (and  perhaps  bias  errors  in  drag  coefficient).  For  the  expanding  19-21 
linear -cubic  smoothing  the  bias  is  given  for  125  km  apogee  in  Figure  5.  Assum¬ 
ing  this  bias  profile,  the  resulting  pressure  error  at  any  point  P.  can  be  evaluated 
as 
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where  Ap.  is  the  signed  (±)  bins  in  density  at  the  ith  point.  The  percent  error 
1  ^i 

in  pressure  — —  resulting  from  the  evaluation  of  &)uation  22  is  presented  in 
i  A 

Table  10.  See  Table  4  for  pertinent  yalues  of  (100)* 
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TABLE  19 

Percent  Error  in  Pressure 

Resulting  from  Bias  Error  in  Density--125  Km  Apogee 
ALTITUDE  ~  KM 
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4.1.4  Constant  Percent  Bias  in  Density 

It  is  interesting  to  observe  the  effect  that  a  constant  percent  bias  in 

Api 

density  has  on  pressure  accuracy.  Assume -  =  K  a  constant.  Then  from 

Pi 

Equation  19 
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Hence  a  constant  percent  bias  in  density  (or  for  that,  matter  in  drag  coefficient) 
produces  the  same  percent  bias  in  pressure. 

4.  2  Summary  of  Pressure  Errors 

Pressure  errors  result  from  the  error  in  the  initial  estimate  of  temper¬ 
ature,  past  and  present  noise  errors  in  density,  and  past  and  present  bias  errors 
in  density,  r  igure  22  shows  the  estimated  magnitude  of  the  resulting  pressure 
errors  for  a  125  km  apogee  using  the  expanding  19-21  lin-cub  smoothing  on  density. 
A  10%  error  in  the  initial  estimate  of  temperature  is  shown  to  diminish  to  les3 
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than  a  1%  error  in  pressure  after  7  km  of  flight.  The  preosure  error  from  the 
noise  in  density  is  less  than  2%  above  90  km  and  less  than  1%  below  90  km.  The 
pressure  bias  introduced  by  the  smoothing  induced  bias  in  density  is  2~3%  above 
82  km  and  decreases  to  less  than  1%  below  78  km.  A  constant  percent  bias  in 
density  produces  the  same  constant  percent  bias  in  pressure 


Altitude  Km 

Figure  22:  Preeeure  Error*  for  Sphere  with  125  Km  Apogee 


5.  TEMPERATURE 

Temperature  is  computed  from  the  pressure  and  density  using  the 
gds  ^aw*  Temperature  errors  result  from  errors  in  density  and  pressure.  In 
addition,  above  90  km  temperature  errors  may  also  result  from  a  change  in 
the  molecular  weight  of  air. 


5. 1  Temperature  Error  Equation 

The  temperature  increment  resulting  from  error  increments  in  pressure 
and  density  is  approximated  by 


dT  = 


5T 
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dp  + dP 
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(24) 


dT  dP  dp 
T  =  P'  '  p 
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The  effect  of  density  and  pressure  errors  on  temperature  depends  upon  the  type 
of  errors  in  density  and  pressure.  Noise  and  bias  errors  occur  in  density.  The 
effect  of  each  of  these  type  errors  on  pressure  accuracy  has  been  described  in 
section  4.  By  combining  each  type  density  error  with  the  resultant  error  it 
produces  in  pressure,  and  substituting  into  Equation  24.  a  temperature  error  profile 
is  generated.  In  addition  to  density  errors  that  produce  pressure  errors,  the 
error  in  the  initial  temperature  also  results  in  pressure  and  temperature  errors. 

5.  L 1  Temperature  Errors  Resulting  from  Error  in  the  Initial  Preseure 

The  error  in  the  initial  estimate  of  temperature  produces  a  pressure  error 
and  aleo  a  density  error  since  the  Mach  and  Reynolds  number-  and  thus  depend 
on  temperature.  The  combined  effect  of  this  pressure  and  density  error  is  given  in 
Table  11. 


TABLE  11 

Percent  Error  in  Temperature 
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ALTITUDE  ~KM 


100 

95 

90 

85 

<85 

Fnm  Table  8 

From  Table  2 

“f-  (100) 

6.0 

1.4 

.6 

.  1 

Lesa  than  3.  1  % 

-y-(ioo) 

-4.0 

8 

-.3 

1 

Less  than  0.  1% 

*^■(100) 

10.0 

2.2 

.9 

.2 

Leas  than  0.  2% 

5.1.  2  Temperature  Errors  Resulting  from  Noise  E'.  in  Density 

Assuming  the  ith  pressure  and  ith  density  errors  to  be  independent, 
(which  is  not  completely  correct  but  they  are  essentially  independent  since  the 
ith  pressure  e.'ror,  after  a  few  data  points,  depends  primarily  on  the  previous 
density  errors  and  only  to  a  small  degree  on  the  ith  density  error),  the  percent 
error  variance  in  temperature  is 


(25) 


The  evaluation  of  Equation  25  using  the  noise  error  density  profile  for 
the  19-21  smoothing  and  the  pressure  error  profile  it  produces  is  presented  as 
Table  12  for  a  sphere  apogee  of  125  Km. 
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TABLE  12 


Percfciu  Error  in  Temperature 

Reoultin.  from  Noire  Enron  in  Density  and  Pressure — 12S  Km  Apogee 

ALTTrVPE  -  KM 


100 

94 

90 

04 

to 

7$ 

70 

*4 

00 

ss 

40 

44 

40 

44 

M 

* 

7rmn  Tabu  4 

« 

Prom  TabU  9 

op 

—  (100) 

0.0 

U 

1.4 

1.0 

0.7 

t.» 

1.4 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

~—{I0OJ 

1.9 

1.1 

0.0 

0.4 

0.  1 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.1 

4.  ? 

0.1 

«T 

— poo) 

0.1 

).  7 

1.7 

1.  1 

0.7 

0.S 

1.  5 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

5. 1,  3  Temperature  Errors  Resulting  from  Bias  Errors  in  Density 

The  bias  errors  in  density  introduced  by  the  19-21  smoothing  and  the 
resulting  pressure  errors  are  given  in  Table  10.  A  substitution  of  these  values 
into  Equation  24  yields  the  resulting  bias  in  temperature.  (See  Table  13). 


TABLE  13 

Percent  Error  in  Temperature 

Resulting  from  Bias  Errors  in  Density  and  Pressure--125  Km  Apogee 


ALTITUDE  -Kl4 


100 

94 

90 

05 

so 

75 

70 

65 

ti 

55 

50 

45 

40 

55 

30 

■"*•(100) 

«*.0 

♦  3.S 

♦  4.0 

♦  3.0 

.  1.  5 

.1.0 

.1.0 

0 

0 

0 

0 

0 

0 

0 

0 

ap  (100) 

p 

,2.0 

♦  2.4 

♦  2.3 

♦  3.1 

♦  1.5 

.0.2 

-0.4 

.  0.2 

.0.  1 

.0.  1 

0 

0 

0 

0 

0 

—  (100) 

♦  0.0 

-  1. 1 

-1.7 

+  0.  1 

3.0 

♦o.e 

_ 

♦0.6 

-  0.2 

-0.1 

•0.  1 

0 

0 

0 

0 

0 

From  Tabla  1° 


5. 1.  4  Constant  Percent  Bias  in  Density 

It  was  shown  in  section  4  that  a  constant  percent  bias  in  density  yielded 
the  same  percent  bias  in  pressure.  Hence,  from  Equation  24  it  follows  directly 
that  no  temperature  error  is  introduced  by  a  constant  percent  eii-or  in  density. 

5.  2  Summary  of  Temperature  Errors 

Temperature  errors  ?rC  the  result  of  errors  in  density  and  pressure 
and,  above  90  km,  may  also  be  the  result  of  a  change  in  the  molecular  weight 
of  air.  Figure  23  shows  the  estimated  temperature  errors  resulting  from  density 
and  pressure  errors  for  a  sphere  with  apogee  125  km.  A  10%  error  in  the  initial 
estimate  of  temperature  is  shown  to  diminish  to  less  than  a  1%  error  in  temperature 
after  10  km  of  flight.  The  noise  error  in  temperature  resulting  from  the  noise 
in  density  is  a  maximum  of  3.  7%  at  95  km  (approximately  8  deg  ),  decreasing 
to  less  than  2%  (approximately  4  deg  )  by  90  km  and  remaining  at  2%  or  lees 


s0 


for  the  remainder  of  the  flight.  The  temperature  bias  produced  by  the  19-21 
smoothing  varieo  between  ±  ?.%  (approximately  4  deg  )  from  100  to  75  km  and  is 
less  than  1/2%  (approximately  1  deg  )  below  75  km.  A  constant  percent  bias  in 
density  produces  no  error  in  temperature. 


o 


Figura  23:  Temperature  Error*  for  iphere  with  125  Km  apogee. 

6.  DESCRIPTION  OF  HIGH  ALTITUDE  ROBIN  PROGRAM 

This  section  summarizes  the  mechanics  of  the  computer  program, 
wirtten  to  implement  the  techniques  derived  in  sections  2-5.  The  program  title 
is,  "May,  1970  AFCRL  High  Altitude  ROBIN  Program.  "  A  users  manual  (Ref.  12) 
which  describes  the  program  and  describes  how  to  interpret  the  output  data 
is  available  through  the  Vertical  Sounding  Techniques  Branch  (CRER)  of  the 
Air  Force  Cambridge  Research  Laboratories.  The  basic  flow  of  the  program, 
order  of  computation,  iteration  procedures,  etc.,  are  akin  to  those  explained’ 
in  Reference  1. 
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6.1 


Editing 


Input  to  the  ROBIN  Program  is  time  and  X,  Y,  Z  position  coordinates 
which  have  been  obtained  by  a  transformation  of  the  radar  observed  slant  range, 
elevation  angle  and  azimuth  angle.  Prior  to  performing  any  computation  on  the 
position  coordinates  an  editing  procedure  is  inaugurated  to  remove  any  spurious 
data  that  may  have  resulted  from  a  radar  malfunction.  The  editing  begins  whan 
the  ascent  part  of  the  rocket  trajectory  exceeds  150,  000  ft.  The  procedure 
consists  of  comparing  the  first  difference  between  consecutive  position  coordinates 
with  a  pre-established  tolerance.  If  the  difference  exceeds  the  tolerance,  the 
second  of  the  two  points  is  considered  bad.  The  succeeding  point  (third  point) 
is  then  differenced  from  the  first  point  and  compared  to  twice  the  tolerance. 

If  the  third  point  is  accepted,  a  linear  interpolation  between  the  first  and  third 
points  replace?  the  bad  second  point.  If  the  third  point  fails  the  tolerance,  the 
procedure  is  repeated  using  the  fourth,  fifth,  etc.  data  point  until  either 
another  good  point  is  found  or  20  consecutive  bad  points  (2  seconds  of  data)  have 
been  determined.  If  20  consecutive  bad  points  are  found,  or  alternately,  if 
2  seconds  of  missing  data  is  found  the  program  is  restarted.  The  tolerance 
used  for  the  X  and  Y  component?  is  150  ft.  This  implies  that  the  maximum 
allowable  velocity  obtained  by  finite  difference  between  two  consecutive  0. 1 
second  data  points  is  1500  ft/sec.  Since  vertical  velocities  often  exceed  2500 
ft/sec  at  high  altitudes,  two  tolerances  are  used  for  testing  the  Z  coordinate. 

The  Z  tolerance  below  200,  000  ft  is  200  ft,  and  above  200,  000  ft  is  360  ft. 

The  editing  procedure,  though  not  overly  sophisticated,  performs  very 
capably  with  typical  ROBIN  flights.  Experience  has  shown  that  the  radar  track 
requires  considerable  editing  on  the  upleg  section  of  the  rocket  trajectory,  and 
again  at  the  very  low  altitudes,  below  30  km,  of  the  balloon  trajectory.  The 
usable  segment  of  data,  between  30  km  and  125  km,  for  computing  meteorological 
parameters  requires  minimal,  and  sometimes  no  editing. 

6.2  Smoothing 

The  program  is  designed  to  accept  FPS-16  Radar  Data  at  0. 1  second 
increments  after  it  has  been  transformed  to  rectangular  (X,  Y,  Z)  coordinates. 

The  program  averages  5  of  the  0.1  second  data  points  to  arrive  at  1/2  second 
data  points  which  are  statistically  independent.  These  1/2  second  spaced 
rectangular  position  coordinates  become  the  input  data  for  the  smoothing 
procedure. 

6.2.1  Density 

Above  60  km  the  vertical  velocity  and  acceleration  required  for  input 
to  the  density  equation  is  computed  by  19-21  linear  cubic  smoothing.  That  is 
19  vertical  position  coordinates  (1/2  second  spaced)  are  fitted  by  a  linear  poly¬ 
nomial  and  the  slope  taken  as  the  vertical  velocity.  A  slide  of  two  data  points 
after  each  fit  results  in  a  series  of  velocities  spaced  one  second  apart  in  time. 

21  of  these  velocities  are  fitted  by  a  cubic  polynomial  and  the  derivative  at  the 


•  _ _ _ 

midpoint  of  the  interval  taken  as  the  acceleration,  A  slide  of  one  velocity  point 
after  each  fit  of  velocity  results  in  acceleration  time  spacing  of  one  second. 

The  horizontal  velocities  and  accelerations  required  for  the  density  equation 
art  computed  by  51-35  cubic -cubic  smoothing.  Below  60  km  the  velocity 
smoothing  interval  is  expanded  in  such  a  way  so  as  to  maintain  an  approximate 
2 %  error  in  density.  This  is  accomplished  by  expanding  the  smoothing  interval 
by  4  points  each  50  seconds.  That  is,  at  time  60  km  +  50  seconds,  23-21  lin-cub 
smoothing  is  initiated  ;  at  time  60  km  4-  100  seconds,  27-21  smoothing,  etc.  , 
until  51-21  smoothing  is  reached.  At  this  point  (approximately  35  km)  the  51-21 
linear -cubic  smoothing  is  retained  for  the  remainder  of  the  flight. 

6.  2.  2  Winds 

All  three  rectangular  coordinates,  vertical  as  well  as  horizontal,  used 
in  the  wind  equation  are  computed  by  51-35  cubic -cubic  smoothing. 

6.  3  Calculation  of  Meteorological  Parameters 

The  derivation  of  the  equations  used  to  calculate  wind,  density,  temp- 
erature„and  pressure  using  the  passive  sphere  technique  is  well  substantiated 
in  Reference  1.  All  forces  acting  upon  the  sphere  are  included  in  the  equation 
with  the  exception  of  the  forces  acting  upon  the  apparent  mass  of  the  sphere.  It 
can  easily  be  shown  that  the  apparent  mass  of  the  sphere  can  be  neglected  with¬ 
out  affecting  measurement  accuracy.  A  flow  diagram  of  the  procedures  used  to 
calculate  the  meteorological  parameters  as  excerpted  from  Reference  1  is  given 
in  Figure  24. 


Figaro  24:  Procedural  Diagram  (or  Calculation  o(  Meteorological  Parameter* 
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6.4 


Determination  of  Density  Validity 


The  validity  of  density  data  is  determined  by  two  checks:  the  time  of 
fall  check  and  the  density  gradient  check  (X  check). 

6.4.1  Time  of  Fall  Check 

The  time  of  fall  check  is  the  primary  test  for  sphere  collapse.  It 

is  used  between  80  km  and  30  km.  To  use  the  time  of  fall  test  above  55  km 

sphere  apogee  must  be  first  calculated.  Having  determined  the  sphere  apogee, 

the  actual  time  of  fall  of  the  sphere  between  80  and  70  km,  between  70  and  t»0 

km,  60  and  55  km,  etc. ,  is  compared  with  that  expected  for  a  sphere  falling 

from  the  determined  apogea,  through  an  atmosphere  of  0.  8  p  ,  ,  and  1.  2  p 

oZ  be 

If  the  actual  time  of  fall  is  not  within  the  range  of  that  expected,  sphere  collapse 
is  indicated.  The  printout  reads  "time  of  fall  test  between  xx  and  xx  km 
indicates  sphere  collapse.  " 

6.  4.  2  Density  Gradient  Check  (Xcheck) 

The  density  gradient  check  is  used  below  60  km  and  overlaps  the  time 
of  fall  test  between  60  and  30  km.  It  is  the  secondary  sphere  collapse  check. 

The  density  gradient(_L  ^L)  is  estimated  below  60  km  by  calculating  X  as  defined 
by  Equation  8.  If  X  defea  not  fall  within  the  range  of 


00014  -  /(.  00005)2  +(  <  X  <.  00014  +J(.00005)2+(- 


08484 y 
Az  / 


sphere  collapse  is  indicated  and  the  printout  reads  "density  gradient  exceeds 
nominal  value. "  It  is  possible  that  a  large  vertical  wind  may  at  times,  trip  the 
X  check.  Thus,  the  X  check  should  be  interpreted  in  conjunction  with  the  time 

of  fall  test. 


6.  5  Measurement  Accuracy 

Included  in  the  calculation  of  each  parameter  (density,  wind,  tempera¬ 
ture,  and  pressure)  is  a  calculation  of  the  estimated  noise  error  in  each  parameter. 
The  estimate  is  made  using  the  formulas:  Density  noise  error  terms  of  Equation 
6;  Winds,  Equation  16;  Pressure,  Equation  21?  Temperature,  Equation  25.  The 
position  error  is  determined  by  assuming  r.m.  s.  errors  of  0,1  mil  in  azimuth 
and  elevation  and  4.  5  yards  in  slant  range.  The  other  types  error  present 
(bias,  error,  etc.  )  are  not  included  in  the  computer  printout.  For  the 
magnitude  of  these  errors  "the  reader"  is  referred  to  Figures  5,  16,  and  23, 
and  Table  1. 
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A  frequency  response  table  for  density  and  wind  is  printed  for  various  escape 
altitudes  at  the  end  of  the  program  (Figure  25/.  This  serves  as  a  means  of 
determining  the  amount  of  detail  observable  in  the  wind  and  density  output. 

6.6  Program  Options  and  Output 

The  program  contains  options  for  the  type  of  output  desired  by  the 
user.  Option  1:  the  program  allows  for  the  output  to  be  piinted  in  either  of 
two  forms;  a)  data  can  be  printed  at  constant  time  increments  of  one  second, 
or  b)  the  data  can  be  printed  to  include  constant  altitude  of  200  meters.  The 
constant  altitude  printout  is  determined  by  linear  interpolation  of  the  constant 
time  data.  Figuref  26  and  27  show  a  sample  of  each  type  output.  Option  2: 
the  program  allows  for  a  plot  of  the  output  data  on  a  Cal-comp  plotter.  The 
temperature,  density  ratio,  pressure  ratio,  and  wind  components  can  be 
plotted  as  a  function  of  altitude.  Figure  28  shows  the  output  data  that  serves 
as  the  input  to  the  plotting  program.  Sample  plots  from  the  Cal-comp  plotter 
can  be  observed  in  Figures  29,  30,  32,  33,  34,  35,  39,  and  40. 

7  OPTIMUM  WIND  AND  DENSITY  REDUCTION  OF  THREE 

OPERATIONAL  FLIGHTS 

A  series  of  three  Viper  Dart  flights  were  launched  at  Eglin  Air  Force 
Base  on  18  February  1968  at  18:00,  19:00  and  20:00  zulu.  Each  flight  was 
tracked  by  two  FPS-16  radars.  The  flights  are  identified  as  Viper  Dart  11,  12 
and  13  and  the  radars  as  radar  23  and  radar  27.  The  flights  were  reduced  using 
the  May  1970  High  Altitude  ROBIN  Program. 

7.1  Density  Accuracy 

Figures  29  and  30  show  the  density  ratio  (  — ^ —  )  for  each  track  of 

p"62" 

each  flight.  For  each  of  the  flights,  there  is  excellent  agreement  between  radar 
23  and  radar  27.  Figure  31  shows  the  percent  noise  error  in  density  obtained 
from  the  printout  of  the  computer  program*.  If  the  lcr  error  estimate  is  valid, 
approximately  2/3  of  the  actual  errors  (density  differences  from  two  radar  tracks 
should  be  within  the  ltr  estimate.  Figure  31  appears  to  substantiate  this  contention, 
especially  below  95  km.  It  is  interesting,  however  to  observe  the  variation  in 
density  between  flights  one  hour  apart  in  time,  particularly  in  the  62-54  km 
region.  The  cause  of  these  differences  is  not  known.  Possible  causes  are: 

a)  An  actual  fluctuation  in  density  between  the  launch  times  of  the 
flights. 


uMrr°rS  err°r  due  t0  the  initial  «uess  of  temperature  are 

essentially  the  same  in  both  radar  tracks.  As  a  result.differencing  densities 

pressures,  temperatures  or  winds  does  not  determine  the  magnitude  of  these’ 
errors. 
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Figure  25:  Program  Output  -  Frequency  Response  of  Density  and  Winl 
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Figure  26:  Program  Output  -  One  Second  Data;Option  la 
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Figure  27:  Program  Output  -  Constant  Altitude  (interpolation);  Option  lb 
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Figure  28:  Computer  Output  of  Data  for  Plotting 
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■VIPER  DART  13  RADAR  27 
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Figure  29:  Density  Ratio  from  Viper  Darts  11,  12,  and  13;  Radar  27 


b)  Spacial  variations  in  dsnsity  between  the  paths  of  three 
spheres. 

c)  Inaccuracies  of  the  drag  coefficients  being  experienced  at 
different  altitudes  for  the  three  spheres. 

d)  A  change  in  the  vertical  motions  of  the  atmosphere. 
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These  discrepancies  in  density,  one  hour  apart  in  time,  are  not  due 
to  the  inability  of  the  radar  to  accurately  track  the  spheres.  Comparison  of 
densities  from  the  two  tracks  of  the  same  sphere  clearly  rule  out  this  cause. 
Nor  are  the  density  discrepancies  thought  to  be  a  result  of  a  sphere  collapse 
or  elongation.  Both  sphere  collapse  checks  indicate  the  spheres  were  still 
inflated  to  at  least  40  km. 
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Figure  30:  Density  Ratio  from  Viper  Darts  11,  12,  and  13:  Radar  23 
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Figure  31:  Comparison  of  Estimated  to  Actual  Density  Errors: 
Viper  Darts  11,  12,  and  13 


7.  2  Wind  Accuracy 


Figures  32-35  are  plots  of  the  W  (East  Wind)  and  W  (North  Wind)wind 
components  for  each  of  the  three  flights.  XAbove  85  km  Viper  ^barts  12  and  13  show 
large  differences  between  winds  from  radar  23  and  radar  27.  To  determine  if 
the  differences  were  within  the  bounds  of  that  estimate,  plots  were  made  compar¬ 
ing  the  actual  noise  errors  to  the  estimated  l<r  errors.  Figures  36  and  37  show 
the  noise  error  in  winds  obtained  by  differencing  the  winds  from  the  two  tracks 
compared  to  that  estimated  by  the  computer  program  output.  Both  the  W  and  W 
component  errors  from  Viper  Dart  11  show  the  expected  agreement.  Viper  Darts  ^ 
12  and  13  show  good  agreement  in  general  only  below  85  km.  The  cause  of  this 
decrease  in  agreement,  which  is  beyond  what  is  predicted  from  the  total  error 
plots  for  51-35  cubic -cubic  smoothing  (Figure  16)  has  been  investigated.  The 
following  results  have  been  obtained. 
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The  tracks  of  radar  23  from  both  Viper  Dart  12  and  Viper  Dart  13  flights 
showed  large  oscillation  in  horizontal  position  which  did  not  appear  in  the  tracks 
of  radar  27  from  the  same  flights  (see  Figure  38).  Previous  experience  with 
FPS-16  tracking  data,  indicates  the  cause  of  the  oscillation  in  the  track  from  radar 
23  is  probably  a  low  servo  band  width  setting.  The  fact  that  radar  27  produces 
a  smooth  non -oscillating  track  indicates  the  oscillations  as  seen  in  the  track  of 
radar  23  are  not  real.  Further  investigation  of  the  effect  of  the  servo  on  tracking 
of  passive  spheres  is  in  order.  In  comparing  the  wind  profile  from  the  flights 
one  hour  apart,  changes  in  the  wind  field  are  seen.  These  changes,  especially 
below  80  km  appear  to  be  real  variations  in  the  wind  field,  since  they  are  observed 
by  both  radars. 
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figure  33:  East  Wind  from  Vip* 
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Figure  35:  North  Wind  from  Viper  Darts  11,  12  and  13:  Radar  23 
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Figure  %•  Comparison  of  Estimated  to  Actual  East  Wind  Errors: 
Viper  Darts  11,  12,  and  13 


Figure  37:  Comparison  of  Estimated  to  Actual  North  Wind  Errors; 
Viper  Darts  11,  12,  and  13. 
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Figaro  38:  Compariaon  of  Y  positions  from  Rod*.?  27  and  Radar  23;  Vipor  Dart  12 


Pressure  Accuracy 


Figures  39  and  40  show  the  pressures  obtained  by  each  radar  for  the 
three  flights.  The  actual  noise  error  in  pressure  (by  differencing  the  pressure 
from  the  two  tracks  of  the  same  flight)  is  compared  to  the  estimated  lo-  noise 
error  in  pressure  obtained  from  the  High  Altitude  ROBIN  Program  output  in 
Figure  41.  The  agreement  between  the  actual  and  the  estimated  errors  is 
good  below  90  km.  Above  90  km  pressure  errors  are  larger  than  that 
estimated.  This  is  apparently  the  result  of  the  low  servo  band  width  setting 
of  radar  23. 
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Figure  39:  Pressure  Ratio  from  Viper  Darts  11,  12  and  13:  Radar  27 
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Figure  40:  Pressure  Ratio  from  Viper  Darts  U,  12,  and  13:  Radar  23 
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Figure  41:  Comparison  of  Estimated  to  Actual  Pressure  Errors;  Viper 
Viper  Darts  11,  I2,and  13, 


7. 4  Temperature  Accuracy 

The  temperature  profiles  for  flights  Viper  Dart  11  -13  are  presented 
in  Figures  42  and  43.  By  differencing  the  temperatures  from  two  tracks  of 
the  same  flight  the  noise  error  in  temperature  is  computed.  Figure  44  shows 
the  actual  noise  error  in  temperature  for  each  flight  and  the  estimated  l<r  noise 
error.  For  each  of  the  three  flights,  the  agreement  is  very  good  except  above 
90  km.  The  disagreement  above  90  km  results  from  the  somewhat  larger  than 
anticipated  density  and  temperature  errors  above  90km. 


71 


40  45  iO 


AAJTTUOf-W 


Figure  42:  Temperature  from  Viper  Darts  11.  12,  and  13;  Radar  27 
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Estimated  le  noiss  erroi 
from  computer  program. 


Actual  noise  error  b»  1  — 

differencing  temperatures  — 
from  radar  27  and  radar  23 
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Conclusions 


A  series  of  three  Viper  Dart  flights  launched  one  hour  apart  has 
been  reduced,  using  the  May  1970  AFCR.L  ROBIN  Program.  The  density 
noise  error  for  these  flights  was  consistent  with  that  predicted  for  the 
expanding  19-21  linear-cubic  smoothing  over  the  entire  altitude  span(30-100  km). 
Observed  differences  in  density  between  the  flights  one  hour  apart  cannot  be 
attributed  to  smoothing  error,  since  the  two  tracks  of  each  sphere  show 
remarkable  consistency.  Wind  accuracy  for  the  three  flights  wap  commensurate 
with  that  predicted  for  the  51-35  cubic -cubic  in  general  only  below  85  km.  Viper 
Dart  11,  however,  showed  good  agreement  over  the  entire  range  from  30  to  100 
km.  The  decrease  in  accuracy  above  85  km  for  Viper  Dart  12  and  13  is  believed 
to  be  caused  by  a  low  servo  band  width  setting  on  radar  23.  Horizontal  position 
plots  from  radar  23  show  oscillations  not  found  in  the  track  of  radar  27.  These 
oscillations  are  often  observed  in  sphere  tracking  at  a  low  servo  band  width. 
They  do  not  appear  predominate  in  vertical  position  plots  (Z,t) indicating 
perhaps  a  low  servo  band  in  azimuth  for  radar  23.  Further  investigation  of 
this  problem  is  necessary.  As  in  the  case  of  density,  the  one  hour  apart 
flights  often  show  large  changes  in  the  structure  of  the  wind  field  below  80  km, 
sometimes  as  large  as  15  m/sec.  Since  both  radar  tracks  (for  each  sphere) 
observe  the  same  wind  detail  (below  80  km)  this  variation  in  wind  structure  is 
believed  to  be  real.  Temperature  and  pressure  accuracy  was  in  good  agree¬ 
ment  with  that  estimated  to  altitudes  of  90  km.  Above  90  km  the  low  servo  band 
width  setting  of  radar  23  is  believed  to  increase  the  errors  in  excess  of  that 
estimated. 


8.  SUMMARY  OF  DENSITY,  WINDS,  PRESSURE  AND 

TEMPERATURE  ERRORS 


For  the  convenience  of  the  reader  this  section  summarises  the  expected 
errors  in  wind  and  thermodynamic  computations  from  the  ROBIN  System  using 
the  May  1970  AFCRL  ROBIN  Program.  Table  14  shows  the  error  estimates  for 
each  of  the  meteorological  parameters  and  the  source  of  these  errors.  All  errors 
listed  in  Table  14  are  point  errors  in  the  sense  that  they  are  the  estimated  errors 
in  a  data  point  at  the  specified  altitude.  In  the  case  of  pressure  and  temperature, 
which  are  computed  from  an  integration  of  density,  the  point  errors  reflect  the 
effect  of  the  errors  in  previous  data  points.  For  each  of  the  meteorological 
parameters  the  noise  error  is  the  error  in  that  parameter  which  results  from 
the  noise  in  the  radar  coordinates.  The  bias  error  is  that  error  resulting 
from  any  bias  introduced  by  the  smoothing  process.  The  error  in  density 
resulting  from  errors  in  the  Drag  Coefficient  are  difficult  to  assess  due  to  the 
inconsistent  results  obtained  by  different  experimenters.  Those  Cp  errors 
quoted  in  Table  14  are  the  UDRI's  best  estimate,  with  the  information  available, 
and  is  subject  to  revision  if  warrented  upon  the  completion  of  work  in  progress 
at  AEDC.  The  error  in  density  from  a  X  m/eec  vertical  wind  is  obtained  by 
multiplying  the  error  from  a  1  m/sec  vertical  wind  (given  in  Table  14)  by  X. 

To  compute  the  effect  of  Cq  errors  on  pressure  and  temperature  the  direction 
(  +  or  -)  of  the  CD  error  as  well  as  its  magnitude  must  be  known.  Since  this  is 
not  known  the  effect  of  the  CD  errors  on  these  parameters  cannot  be  determined. 
To  compute  the  effect  of  vertical  winds  on  pressure  and  temperature  a  vertical 
wind  profile  vs  altitude  must  be  assumed.  Since  little  is  known  of  the  structure 
of  vertical  winds  in  the  atmosphere  a  representative  profile  is  not  available. 

The  noise  error  in  winds  can  be  thought  of  as  the  error  in  the  mean  wind  profile.. 
An  estimate  of  the  bias  error  in  wind  is  not  given  because  it  depends  upon  the 
wind  field  itself.  For  a  given  wave  structure  of  the  wind  field  the  distortion  of 
any  frequency  resulting  from  the  smoothing  process  is  given  by  the  frequency 
response.  Table  15  shov/s  the  frequency  response  for  both  wind  and  density  data. 
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SUMMARY  OF  DEfSITY,  PRESSURE.  TEMPERATURE.  AND  WIND  ERROLS 


Altitude  (Km) 

MM 

45 

40 

as 

•0 

75 

70 

4$ 

to 

5$ 

50 

40 

10 

Noise  Error 

• 

1.5 

L  5 

l 

0.7 

0.0 

L  5 

2 

2 

2 

2 

2 

2 

s 

Biss  Error 

+2 

+1.5 

*4 

+3 

-LS 

•1 

-1 

teas  than  1/2% 

o 

<*? 

H  • 

Error  doe  to 

10%  error  lo 

4 

0.1 

less  than  1/2% 

.  u 

T 

Eft 

U 

Error  (romCp 

19 

• 

« 

3 

1 

2 

6 

4 

3 

2 

2 

2 

2 

a 

Error  from  T 
1  m/sec  |L.± 

T«i 

.01 

.01 

.01 

.01 

.02 

.03 

.10 

.3 

.7 

1 

4 

4  ^ 

vertical  « 

.01 

.01 

.01 

.01 

wind  ^£2. 

.01 

.2 

.4 

.7 

L6 

L  5 

2.0 

5.0 

Xo.o 

*  10 

•  2 

.2 

.2 

.2 

.2 

.4 

.5 

.0 

L0 

L  5 

2.0 

5.0 

10.0 

Noise  Error 

1 

L2 

0.7 

less  than  1/2% 

s 

Bias  Error 

+2 

♦2.4 

+2.3 

♦3.  1 

+L5 

lass  than  1/2% 

o 

PS 

Error  doe  to 

W  o 

10%  error  In 

6 

L4 

0.6 

less  than  1/2% 

S2 

Si 


Error  due  to 

Op  error  and  Cannot !»  determined 
vertical  wind* 


H  2 
<  » 
2  o. 
w  - 

04 

S 

B 


Noiee  Error 

8.1 

1.7 

L  7 

LI 

0.1  0. 8 

1.5 

2  2  2 

2 

2 

2 

Biae  Error 

0 

•1  1 

-1.7 

+  0.  1 

+3  +0.8 

40.6 

lass  than  1/2% 

Error  due  to 

c 

e 

10%  error  in 

10 

2.2 

0.4 

less  than  1/2% 

Error  duo  to 
Cp  e  ror  and 
vertical  wind* 


CD  e  ror  and  Cannot  be  determined 


w  — 
q  55  O 

a28 

*sl 


Noiee  Error 


13 


lets  than  1  m/oec 


77 


Tub!*  15 


Frequency  Rtipoui  of  Wind  uni  Oiutty  i 


125  Km  Apogoo 


«»tic  cf  affutuce  of  smoothed  density  wave  tc  <mirit(  of  cricimi 

M«V!  AS  A  Ft'iCUCW  CF  AllHbDE  AKO  WAVELENGTH. 

(AIT IllOt  AND  WAVELENOTl-  F£A$U«E:  IK  KIICFET£M$.|  ■ 

AIT  mot 

x  mo  x  «o  x  tc  x  ?c  >  e«  >  »c  x  40  x  so  x 

xxxxxxxxxxxxvxxxxxxxxxxxxxixmxxxxxxxxxxxxxxxxxxxxxxxxxx 
■  xxxxx  xxxx 
I  1  X  C.01  X  C.C1  X  C.C1  X  C.C1  X  C.Cl  X  c .22  X  C.4i  X  c.t»  X 

rxxxxxxxxx 

2  X  0.01  X  C.01  X  C.C1  X  t.CA  X  C.22  X  0.12  X  C.tC  X  C.90  X 

xxxxxxxxx 
5  X  0.C1  <  C.C2  X  C.OJ  X  (.47  X  C.d4  X  C.40  X  0.46  X  (.49  X 

IX  X  X  X  X  X.X  X  X 

i  10  X  0.67  X  C.S1  X  0.09  X  (.72  X  C.47  X  C.98  X  C.99  X  C.99  X 

t  .X  X  X  X  X  X  X  X  X 

I  20  X  0.94  x  C.99  X  0.96  X  C.SS  X  C.99  X  C.99  X  0.99  X  C.99  X 


RAT  1C  CF  AOFLlTLOt  OF  SMUCIhEC  SIKUSQICAl  MIKQ  TO  AFFUTUCE  Of 
ORIGINAL  W INC  AS  A  FUNCTION  OF  A(T1TU0C  AKC  WAVEIENCTM. 

(ALT1TUCE  AflC  WAVELENGTH  FtASUREO  IK  KIKFETERS.l 

ALT1H0E 

X  ICO  X  40  X  tO  ?  7?  X  tO  X  5C  X  4C  X  30  X 

XXXXXXXXXXXXXXX XXX XX XXXXX XXXXXXXXX) XX XXXXXXXXXXXXXXXXXXXX 
RXXXXXXXXX 
A  IX  0.01  X  C.C1  X  C.01  x  C.d  X  C.01  X  C.C2  X  0.97  X  C.96  X 

V  X  X  X  X  X  X  X  X  X 

E  2  x  c.01  X  C.01  X  C.01  X  C.C3  X  C.1C  X  C.40  X  C.76  X  C.99  X 

(XXXXXXXXX 
E  9  X  C.52  x  C.01  X  c.ce  X  C.12  X  C.7C  X  0.92  X  C.99  X  C.99  X 
NXXXXXXXXX 
9  10  X  0.14  X  C.17  X  C.24  X  C.29  X  C.97  X  C.99  X  0.99  X  C.99  X 

TXXXXXXXXX 
h  20  x  C.56  X  C.63  X  c.tt  X  C.TC  X  C.9e  X  0.99  X  C.99  X  (.99  X 
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APPENDIX  A 


distribution  is  given  by 


Denaity  Error  Due  to  Errors  in  V. 

The  increment  of  error  in  denaity  for  horizontal  velocity  errors  of  e. 

and  e.  and  horizontal  wind  errors  of  e  and  e  is  given  by  - 
y  w  w  ®  1 

x  y 


e  =  -&-C  =  +  •Si-c  +  . 

p  av  v  aJc  1  a*  y  aw  aw  wy 

x  y 


the  percent  error  increment  is: 


X  -  W 


Y  -  W 


.2  x 


-  &»  + 
r2  y 


X  -  w 

- r-^-C  + 

V2  Wx 


Y  -  W 


2  w 
V  y 


Since  the  errors  inX  and  Y  are  much  smaller  than  the  errors  in  W  and  W  , 

x  y 

the  primary  contribution  to  Eq.  1A  comes  from  the  last  two  terms.  Hence, 
Eq.  1A  can  be  approxima.  ^  uy 


X  -  W 

x 

- = - e  + 

z  w 
V  x 


Y  -  W 


Assuming  normally  distributed  and  independent  errors  in  and  W  the 
percent  error  variance  in  density  becomes 


<r  .  2  _  X  -  W  2 

A  \  ^  V 


(-^)  ■  L 


(T  + 

W 

X 


Y  -  W  2 

[— 


A-2 


The  assumptions  e  .  <<  e  and  e  .  <<  e  are  not  critical  m  the  derivation 
x  w  y  w 

x  y 

of  an  error  equation.  Even  if  the  assumptions  introduced  a  10%  error  in  the 
error  equation,  the  results  would  be  satisfactory  for  our  purpose.  The  impor¬ 
tant  criterion  is  to  measure  the  meteorological  parameters  as  accurately  as 
possible  and  to  have  an  estimate  of  their  errors.  If  the  estimate  is  off  by  10%,  it 
is  of  little  consequence. 

Density  Error  Due  to  Vertical  Velocity  and  Acceleration  Errors 

The  error  increment  in  density  resulting  from  vertical  velocity  and 

acceleration  errors  of  e.  and  e.  .is  given  by 

z  z 


e 

P 


JKL_e 

••  V 
az  z 


(2A) 


Velocity  and  acceleration  are  computed  by  smoothing  vertical  position  coordinates. 

The  smoothed  velocity  and  acceleration  contain  noise  error  due  to  radar  noise 

and  bais  error  introduced  by  the  smoothing  process.  Interpreting  Eq.  2A  in 

terms  of  noise  error  in  velocityte.l  and  acc deration le..^  the  percent  noise 

z  z 

error  variation  in  density  is  given  by 


ep  V2  +  (Z-w/  1 

“  "  2  •  e,z  +  ez 

P  V  (Z  -  W  )  z  Z  -  g 

z 


(3A) 


By  assuming  independent  and  normally  distributed  errors  in  velocity  and 

•  m 

acceleration  and  observing  that  W  <<Z  and  V*Z  Equation  3A  simplifies  to 

z 


A-3 


The  squared  bias  error  in  density  introduced  by  bias  velocity  and  acceleration 


errors  is  given  by 


where  Ap  represent  the  bias  error  in  the  parameter  p,  and  Az  and  Az  the  bias 
•  •  • 

in  Z  and  Z.  The  total  percent  error  in  density  is  defined  as 
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total 


or 


<r.. 
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e  e 

Z  -  g 
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a  •• 

Az 

•  e 

Z  -  g 


Density  Error  Due  to  Vertical  Winds 

The  percent  variation  in  density  resulting  from  a  vertical  wind  error  of 
A«z  is,  to  a  first  order  approximation 

2  2 
AO  V  +  (Z  -  W  ) 

=  —  —  Aw  .  (4A) 

p  V  (Z  -  W  ) 

z 

•  e 

Since  <  <  Z  and  Vsf  Z,  the  squared  percent  error  in  density  for  a  vertical 

wind  error  of  A  w  is 
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(-/=  (^)2- 


A -4 


‘****+*f >e 


APPENDIX  B 

DERIVATION  OF  THE  NOISE  ERROR  FORMULAS 


a)  Velocity 

a 

Let  X(i)  be  the  radar-observed  data  point  at  time  t  =  i,  X(i)  the  true  value 
of  the  data  point  at  t  =  i  and  c  (i)  the  error  in  the  observed  data  point;  thus 
X  (i)  =  X(i)  +  €  (i). 

If  N  equally  time  spaced  data  points  are  fitted  by  least  squares  to  a 

polynomial  of  degree  k,  then  the  smoothed  velocity  at  the  midpoint  of  the  inter¬ 
im  +  1 

val  ^ tii  point)  is  given  in  terms  of  Legendre  Polynomials  as 


where 


*  N  +  l  1  N  +  l  '  N  +  l  '  N  +  l 


N 

V 

L  P.  (i)  X(i) 


i=l 


At 


Pk2(i) 


i=l 


P  (i)  =  1, 
o 

_  N  +  l 

P1  (l)  =  1  -  2  ’ 


Pk+1 (l)  =  PlPk  ' 


k  (N2  -  k2) 
4(4  k2  -  1) 


k  -  1  ’ 


At  =  time  spacing  between  data  points, 


N  -  Number  of  data  points  used  in  smoothing, 


B  -l 


k  =  Degree  of  smoothing  Polynomial 


P^tt)  =  First  derivative  with  respect  to  i  of  P^fi). 


Substitution  of  X(i)  +  e  (i)  for  X(i)  in  each  A  gives 
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£  Pjtijfxti)  +  e(i)] 


l  Pk(i)  [*»  +  c(i) ] 


.  «L 


N 

At  r  „  2 


l  pi  » 


P  (jl«,  + . . .  +  id - P; (n±l,. 

2  At  9  2  2 


Rearranging  in  termr  of  e  (i)  yields 


N  1  N+l 

P  (i)  P  (—  ) 

-  l  {  t  /  2  +  • 

2  t=i  it  ?  p2(i> 


1  N+l 

.  *  -  [an  *  .»]■ 
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Assuming  that  the  distribution  of  each  c  (i)  is  normal  with  mean  C  and  variance 

c and  that  cov  (e  (i),  e  (i) )  =0  for  i  t  j  the  error  variance  of  the  smoothed 
x 

velocity  is 


<r  f  ,  N+l 

x  (— ) 


N  _  ...  '  .  N+l  . 

V  S  P1  1  P1  (  2  ) 
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Since  ^  Pj  (i)  pk  (i)  =  0  for  j  /  k  the  above  equation  simplifies  to 
i=i 


<r  .  JN+1 . 
x(— ) 


r  '  ,  n+i  ,-i‘ 

f  LP1  (~~2~^j 


r  '  .  N+l,i 
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At 
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It  can  easily  be  shown  that  for  even  k 


'  ,  N+l.  n 

P.  j  •  (  r  )  ~  0 . 
k(cven)  2 


Hence  equation  (IB)  simplifies  to 
_  ■  «• ,  *  _ 2  _  , 
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i=l 


The  terms  in  parenthesis  on  the  right  side  of  the  above  equation  can  be  written 
as  a  function  of  N.  (See  reference  10).  The  expressions  for  the  first  two  terms 


arc 


[p»] 


12 


N 

l 

i=l 


pf(i) 


N  (N  -1) 


(3B) 


B-3 


and 


7(3N2-7)2(N-4)  l 
(N  +  3)  1 


(4B> 


By  substitution  of  3B  and  4B  into  2B,  the  noise  error  in  velocity  is  determined 
for  polynomial  smoothing  of  degrees  1  through  4.  Note  that  since 


P_  )  =  P A  =  0  both  the  velocity  obtained  from  a  second  or 

2  2  4  2 

forth  degree  polynomial,  and  the  noise  error  in  this  velocity  are  identical  to 
that  obtained  from  a  first  or  third  degree  polynomical  respectively. 


Acceleration 

The  noise  error  in  acceleration  for  double  polynomial  smoothing  is 
obtained  by  two  applications  of  the  appropriate  velocity  noise  error  formulas. 

For  accelerations  that  are  obtained  as  the  second  derivative  of  the  least 
squares  fit  of  a  polynomial  of  degree  k  over  N  data  points,  the  noise  error 
formulas  are  derived  as  follows:  The  smoothed  acceleration  at  the  midpoint 
of  the  interval  is  given  by: 
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N+l 
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(5B) 


where  P'^  is  the  second  derivative  of  P^.  By  writing  X  (i)  =  X(i)  +  e (i)  and 
rearranging.  Equation  5B  becomes 
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The  error  Variance  of  the  smoothed  acceleration  is 
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It  can  be  shown  that  for  odd  k 


”  .N+l.  A 

Pk(odd)  2  =  °* 


This  Equation  (6B)  simplifies  to 
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The  first  two  terms  in  parenthesis  when  written  in  terms  of  N  are 


(8B) 
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and 


N  , 
IP4<‘» 
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44100  (3N2-13)2 _ 

2  -  1)(N6-29N4  +  244N2.  576) 


(9B) 


A  substitution  of  Equationt(8B)  and  (9B)  into  (7B)  will  yield  the  noise  error 
formulas  for  acceleration  obtained  as  the  second  derivative  of  a  Quadratic 
and  Qiartic  Polynomial  fit. 
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APPENDIX  C 

BIAS  ERROR  FORMULAS 
Velocity 

Assume  true  position  is  given  by  a  k  th  degree  polynomial 

that  is 

X(t)=a  +a,t+a  ,t2+. . . +a  tk. 

O  I  b  K 

Initially  we  will  derive  the  bias  error  for  velocity  obtained 
as  the  first  derivative  of  a  linear  polynomial  fit  to  the  true  position. 

Velocity  from  Linear  Polynomial  Fit 

Let  N  (odd)  be  the  number  of  position  points  used  in  the  fit. 


N+l 

Let  t  correspond  to  the  time  at  the  midpoint  ( — r-th  point)  of  the  smoothing 
o  2 


system.  If  the  coefficients  of  the  k  th  degree  position  field  are  given  with  re¬ 
spect  to  the  coordinate  system  with  origin  at  tQ  then  we  have  tQ=0.  )  The  rela¬ 
tionship  between  t  and  numbering  of  the  data  points  by  i=l,  2,  . .  .  N  is 
N+l 

t=t  -(——-)  At  +  iAt. 
o  Z 


C-l 


In  terms  of  Legendre  Polynomials  the  velocity  at  the  midpoint  of  the 
smoothing  interval,  for  linear  smoothing  is 


Pjtf)  X(i) 


X{*~-)  =  — 


4t  ?pf 


or  by  substitution  for  x(i) 
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vhere 


X(i)  =  true  position  at  time  t  =  tQ  -  — )  At  +  iAt, 


„  N+l 

Pj(»)  =  i  "  — —  • 


PjU)  =  1. 


Since  t  =  t  -  ( — r — )  At  +  iAt  for  convenience  t  can  be  written  as 
o  c 


t  =  t  +  P.  (i)  At. 
o  l 


C  -2 
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Substituting  (2C)  into  (1C)  yields 
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The  terms  Inside  brackets  in  Equation  ( 3C)  are  evaluated  by  use  of  the 
following  identities  characteristic  of  Legendre  Polynomials: 
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iil  i=l 
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t  pi3,i)=0 
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The  first  few  terms  of  Equation  (3C)  in  brackets  then  become 
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Substituting  these  terms  into  (3C)  gives 


*  (¥>  Y2Vot3VoH  4^^ry- ■  • 


The  true  velocity  at  the  midpoint  is 


X(*0)=V2Vo+3Vo2+4a4*o3+- 


Thus,  the  bias  error  in  velocity  (defined  as  the  difference  between  the  true 
and  fitted  velocities)  is  given  for  linear  velocity  smoothing  by 


■ »» it2  (■ +  4a4  *o M 


2(iV1-)4--- 


(4C) 
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If  the  coefficients  *q,  a^. . .  are  those  associated  with  the  coordinate  system 
with  origin  at  i  (L  e. ,  t  =0)  (4C)  simplifies  to 


AX 


.  2  /  3N  -  7  \ 

X  =  a^At  y — — - j  +  terms  in  a^,  a^»  ... 


VELOCITY  FROM  CUBIC  POLYNOMIAL  FIT 


The  cubic  velocity  fit  of  a  kth  degree  position  field  is  given  by 
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where 
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Recalling  that  the  first  term  of  Equation  (5C)  was  evaluated  above.  Equation  (5C) 


can  be  rewritten  as 
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The  following  identities  Are  useful  In  evaluating  the  terms  in  brackets. 
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£  PjJiJPjliJ'-O, 
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2 

Using  these  identities  the  terms  in  brackets  of  Equation  (6C)  become 
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Substituting  into  (6C)  yields 

XVT)  *  V*Vo  +  3  Vo  +  +  ... 

The  bias  error  in  velocity  for  the  cubic  fit  is 
Ax  =  terms  in  ag,  a^. 

If  the  coefficients  a  ,  a  ,  .  .  .  are  associated  with  the  coordinate  system 

o  1 

with  origin  tQ  =  0  the  bias  error  formula  simplifies  to 
Ax  =  terms  in  a,.,  a^,  .  .  . 

If  velocity  is  obtained  by  a  second  or  fourth  degree  polynomial  fit,  then 
this  bias  in  velocity  is  the  same  as  that  for  a  linear  or  cubic  fit  respectively. 
This  follows  because  the  slope  of  a  linear  and  quadratic  fit  are  identical  at  the 
midpoint  of  the  interval.  Similarly,  the  slope  of  a  cubic  and  quartic  fit  are 
identical  at  the  midpoint  of  the  interval. 
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Acceleration  trom  second  derivative  of  quadratic  fit 

The  acceleration  obtained  from  second  derivative  of  quadratic  fit  la 

calculated  a  a 
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where 
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By  writing  the  position  as  a  kth  degree  polynomial  and  substituting  into  (7C)  givea 


- 


At 


2^ 


L,  P22(1) 

1=1 


H 

\  )  P. 
O  Lj  < 

i=l 


:(i)  +*ilP2{l)(fo  +  pi(lMt)+' 

i=i 


•+Op 


k^2(1)(to  +  P1<1>  At> }  (8C) 
i=l 


By  using  the  identities 
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equation  (8C)  reduces  to 
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The  bias  in  acceleration  is 


Ak*  =  2a. A  At2  )  +  10a5toAt2(-^~^  )t 


(9C) 


If  the  origin  of  the  coordinate  system  is  t  =  0  then  Equation  (90  simplifies  to 


„  ,  2,  3N  -13  . 

Ax  =  2a  ^  At  ( - — - )  +  terms  in  ag,  ag,  .  .  . 


C-9 


The  derivation  of  the  bias  error  for  the  second  derivative  of  a  quartic  fit 

of  a  kth  degree  polynomial  is  derived  in  a  manner  similar  to  the  above 

derivation.  The  first  term  of  the  bias  error  expression  will  contain  a^  and 

t  .  If  the  origin  of  the  coordinate  system  is  t  =  0  then  the  bias  error  formula 
o  o 

will  consist  of  terms  that  contain  only  the  even  coefficients  a^,  a^,  .  .  . 
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